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Abstract 

In the present lectures the foHowing topics are considered: general properties 
of neutrinos, neutrino mass phenomenology (Dirac and Majorana masses), neutrino 
masses in the simplest extensions of the standard model (including the seesaw mech- 
anism), neutrino oscillations in vacuum, neutrino oscillations in matter (the MSW 
effect) in 2- and 3-flavour schemes, implications of CP, T and CPT symmetries for 
neutrino oscillations, double beta decay, solar neutrino oscillations and the solar neu- 
trino problem, and atmospheric neutrinos. We also give a short overview of the results 
of the accelerator and reactor neutrino experiments and of future projects. Finally, we 
discuss how the available experimental data on neutrino masses and lepton mixing can 
be summarized in the phenomenologically allowed forms of the neutrino mass matrix. 
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1 Introduction 



In June of 1998 a very important event in neutrino physics occured - the Super-Kamiokande 
collaboration reported a strong evidence for neutrino oscillations in their atmospheric neu- 
trino data This was not a thunderstorm in a clear sky: the evidence for oscillating 
neutrinos has been mounting during the last two decades in the solar neutrino experiments 
0, and previous data on atmospheric neutrinos also gave indications in favour of neutrino 
oscillations In addition, the LSND Collaboration has reported an evidence for z/^ v^. 
and z/^ ^ z/e oscillations in their accelerator neutrino experiment However, it was the 
Super-Kamiokande experiment that, for the first time, not only showed with a high statistics 
the deficit of the detected neutrino flux compared to expectations, but also demonstrated 
that this deficit depends on the neutrino pathlength and energy in the way it is expected 
to depend in the case of neutrino oscillations. Since neutrinos are massless in the standard 
model of electroweak interactions, the evidence for neutrino oscillations (and therefore for 
neutrino mass) is the first strong evidence for physics beyond the standard model. 

Neutrino physics is a very active field now, both experimentally and theoretically. The 
Sudbury Neutrino Observatory (SNO) experiment was put into operation in 1999; it should 
add a very important information to our knowledge of solar neutrinos. The SNO experiment 
and the forthcoming Borexino experiment will complement the already existing data of 
Homestake, Gallex, SAGE, Kamiokande and Super-Kamiokande experiments and this will 
hopefully lead to the resolution of the long-standing solar neutrino problem. The first long 
baseline terrestrial neutrino experiment, K2K, has already started taking data, and several 
more experiments - KamLAND, MINOS, MiniBooNE and CERN - Gran Sasso experiments 
are either scheduled to start in a very near future or in advanced stage of planning. There are 
also very wide discussions of possible next generation long baseline experiments with muon 
storage rings ||^. All these experiments are designed to probe a wide range of the neutrino 
mass squared differences and lepton mixing angles, and possibly CP-violating effects in 
neutrino oscillations. They may also be able to test the fascinating possibility of matter 
enhancement of neutrino oscillations - the MSW effect 0. 

On the theoretical side, there have been many analyses of the available experimental data 
and predictions for the forthcoming experiments made within various scenarios for neutrino 
properties, studies of neutrino-related processes in stars and supernovae, phenomenological 
studies of the allowed structures of neutrino mass matrices and developments of particle 
physics models capable of producing the requisite mass matrices. 

Why are we so much interested in neutrinos? Neutrinos play a very important role 
in various branches of subatomic physics as well as in astrophysics and cosmology. The 
smallness of neutrino mass is very likely related to existence of new, yet unexplored mass 
scales in particle physics. These scales are so high that their direct experimental study 
may not ever be possible. Neutrinos can provide us with a very valuable, though indirect, 
information about these mass scales and the new physics related to these scales. In fact. 
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they may even hold a clue to the general problem of the fermion mass generation. 

There are intricate relationships between neutrino physics and nuclear physics. The 
existence of neutrinos was first theoretically suggested to explain apparent energy non- 
conservation and wrong spin-statistics relations in nuclear beta decay; inverse beta decay 
processes were used for the first experimental detections of reactor antineutrinos and solar 
neutrinos. The chiral nature of neutrinos, parity nonconservation and V — A structure of 
weak interactions were also established in nuclear physics experiments 0]. Neutrino-nuclear 
reactions serve as a very clean probe of nuclear and nucleon structure; on the other hand, 
nuclear physics provides us with the knowledge of the reaction cross sections which is very 
important for calculating the solar neutrino fluxes and neutrino detection rates. 

Neutrinos play a very important role in astrophysics and cosmology. They carry away up 
to 99% of the energy released in the type II supernova explosions and therefore dominate the 
supernova energetics. Neutrino reactions play a crucial role in the mechanism of supernova 
explosions. Neutrinos are copiously produced in thermonuclear reactions which occur in 
the stellar interior and in particular in our sun. Solar neutrinos carry information about 
the core of the sun which is unaccessible to direct optical observations. The detection of 
solar neutrinos has confirmed the hypothesis that the sun is powered by thermonuclear 
reactions. At the same time, the sun and supernovae give us a possibility of studying 
neutrino properties over extremely long baselines and probe the neutrino mass differences 
as small as 10~^ eV or even smaller, beyond the reach of the terrestrial neutrino experiments. 

The big bang nucleosynthesis depends sensitively on neutrino interactions and on the 
number of light neutrino species. Neutrinos of the mass m,y ~ 1 eV could constitute the 
so-called hot dark matter and may be important for galaxy formation. Neutrinos may also 
play an important role in baryogenesis: the observed excess of baryons over antibaryons in 
the universe may be related to decays of heavy Majorana neutrinos. 

In the present lecture notes a number of issues pertaining to neutrino physics are consid- 
ered. These include general properties of neutrinos, neutrino mass phenomenology, neutrino 
masses in the simplest extensions of the standard model (including the seesaw mechanism 
of neutrino mass generation), neutrino oscillations in vacuum, neutrino oscillations in mat- 
ter (the MSW effect), implications of CP, T and CPT symmetries for neutrino oscillations, 
solar neutrino oscillations and the solar neutrino problem, and atmospheric neutrinos. We 
also give a short overview of the results of the accelerator and reactor neutrino experiments 
and of future projects. Finally, we discuss how the available experimental data on neutrino 
masses and lepton mixing can be summarized in the phenomenologically allowed forms of 
the neutrino mass matrix. Except for the simplest extensions of the standard models, we do 
not discuss the models of neutrino mass since they were covered in the lectures of K.S. Babu 
at this school (see also |p for recent reviews). Astrophysical implications are limited to 
the solar neutrino problem; other issues of neutrino astrophysics were discussed in the lec- 
tures of T.P. Walker |]T0[, and additional information can be found, e.g. in |Tl|, |T2|, |13[. We 



do not discuss neutrino electromagnetic properties and their implications because of the lack 
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of space. Some material on these issues can be found in [|l^, |T5|, 0. Informative reviews 
and monographs on neutrino physics include (but are not limited to) |T^, |T^, [IS], |T^, ^ . 



2 General properties of neutrinos 

What do we know about neutrinos? Neutrinos are electrically neutral particles of spin 1/2. 
There are at least three species (or flavours) of very light neutrinos, Ue, and z/,-, which 
are left handed, and their antiparticles Ue-, and Ur, which are right handed. Electron- type 
neutrinos and antineutrinos are produced in nuclear (3^ decay. 



A{Z, N) ^ A{Z + 1, iV - 1) + e" + z7e 
A{Z,N) ^ A{Z -l,N + l) + e+ + ue, (1) 

and in particular in the neutron decay process n —>■ p + e~ + Ue- They are also produced 
in muon decays /x^ — > + i^fj.{i^^) + t'e(t'e), and as a subdominant mode, in pion decays 
tt"^ —>■ -\- z/e(t'e) and in some other decays and reactions. The elementary processes 
responsible for the nuclear beta decays or pion decays are actually the quark transitions 
u —>■ d + + i>e and d —>■ u + e~ + Ue- Muon neutrinos and antineutrinos are produced 
in muon decays, pion decays 11'^^^'^ + v^iv^) and some other processes. Neutrinos 
of the third type, tau, are expected to be produced in decays. They have not been 
experimentally detected yet Q, but there is no doubt in their existence. Neutrinos of each 
flavour participate in reactions in which the charged lepton of the corresponding type are 
involved; these reactions are mediated by bosons. Thus, these so-called charged current 
reactions involve the processes ^ + ^a{T^a) where a = e, or r, or related processes. 
Neutrinos can also participate in neutral current reactions mediated by Z^ bosons; these 
are elastic or quasielastic scattering processes and decays Z^ i^ai^a- 

The latter process allowed us to count the number of light neutrino species that have 
the usual electroweak interactions. Indeed, neutrinos from the Z° decays are not detected, 
and therefore the difference between the measured total width of the Z^ boson and the sum 
of its partial widths of decay into quarks and charged leptons, the so-called invisible width, 
Tinv = Ttot — r«s = 498 ± 4.2 MeV, should be due to the decay into pairs. Taking into 
account that the partial width of Z^ decay into one uu pair r,yp = 166.9 MeV one flnds the 
number of the light active neutrino species : 



r 

= = 2.994 ± 0.012 , (2) 

up 

in a very good agreement with the existence of the three neutrino flavours. There are also 
indirect limits on the number of light (m < 1 MeV) neutrino species (including possible 

-'^The DONUT experiment at Fermilab is looking for i^r and at present has several candidate events. 
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electroweak singlet, i.e. "sterile" neutrinos) coming from big bang nucleosynthesis: 

N, < 3.3 , (3) 

though this limit is less reliable than the laboratory one (Q), and probably four neutrino 
species can still be tolerated 




Figure 1: Kurie plots for rrijy = (solid line) and m^y ^ (dashed line) 
What do we know about the neutrino mass? Direct kinematic searches of neutrino mass 



produced only the upper limits |£l 



m^i < 2.5 eV at 95% c.l. (Troitsk) ; < 6 eV at 95% c.l. (Mainz) 
< 170 keV at 90% c.l. (PSI; 7r+ ^ /i+ + z/^) ; 
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< 18.2 MeV at 95% c.l. (ALEPH; ^ 57r + v^) 



(4) 

(5) 
(6) 



Here z/i, z/2 and z/3 are the primary mass components of Ve-, and v-r-, respectively. The limits 
in eqs. (|^) and (|^) come from the comparison of the total energy release with the energy of 
decay products. The limits in eq. @) are obtained from the tritium beta decay experiments 
and are based on the analyses of the so-called Kurie plot. The electron spectrum in the 
allowed /? decay is 



[m,. 



0) 



N,{E,)dE, oc F{Z,E,)^El-mlE,{E^-E,fdE,, 
N,{E,)dE, oc F(Z, E,)JeI - ml E,{Eo - E,)J{Eo - E,f - ml dE, 



(7) 

(m, ^0).(8) 



Here F{Z, Eg) is the well known function which takes into account the interaction be- 
tween the emitted electron and the nucleus in the final state and Eq is the energy re- 
lease. Thus, the shape of the jS spectrum should depend on whether or not neutrinos 
have a mass. It follows from these equations that the plot of the Kurie function K{E(.) = 



4 



[Ne{Ee) /{F(Z, Ee)peEe)Y^'^ veisus Ee should be a straight line when rrii, = but should have 
a different shape close to the endpoint of the spectrum when mj, 7^ (fig. |I|). Analyses of 
the Kurie plots for beta decay of tritium which has a very low energy release give the results 
of eq. (^). However, all the experiments showed some excess of the number of electrons 
near the endpoint of the spectrum rather than a deficiency that is expected if 7^ 0. This 
excess is most likely due to unknown systematic effects. For this reason the conservative 
upper limit on the lightest neutrino mass recommended by the Particle Data Group |21] is 

mi < 15 eV (9) 

rather than the limits in (^). 

Are neutrinos Dirac particles (like quarks or charged leptons) or Majorana particles? In 
other words, can neutrinos be their own antiparticles? It is known experimentally that neu- 
trinos emitted in (3~ decay (which we call electron antineutrinos) cannot be captured in re- 
actions which are caused by electron neutrinos; for example, the reaction i/e+^'^Cl^'^''Ar+e^ 
does not occur, whereas the reaction z/e+^''Cl^'^^Ar+e~ does, and in fact was used in the 
first experiment to detect the solar neutrinos. Does this mean that neutrinos are different 
from antineutrinos, i.e. cannot be Majorana particles? Not necessarily. The reason is that 
the weak interactions that are responsible for neutrino interactions are chiral {V — A): in 
the Cl-Ar reaction mentioned above, only neutrinos of left handed chirality can be detected. 
The particles which we call Ue are right handed, i.e. have the "wrong chirality", and so 
cannot participate in the Cl-Ar reaction (this is called "chiral prohibition"). If neutrinos 
are massive, the chirality (i.e. 1 ± 75) is not a good quantum number, and an electron 
antineutrino which was produced right handed can develop a small left handed component, 
so that the reaction i^e+^^Cl^^''Ar+e~ can become possible. However, the "wrong chi- 
rality" admixture is of the order rrii^/E, i.e. extremely small due to the smallness of the 
neutrino mass compared to typical beta decay energies. Therefore the probability of the 
z/g-l-^^Cl— s>^'^Ar+e~ process should be suppressed by the factor {m^/E)'^, which can be the 
reason for its non-observation. 



Vg ' ^ ' + 




Figure 2: Charged pion decay 



On of the most outstanding manifestations of the chiral prohibition rule are the decays of 
charged pions. They also show that one should carefully discriminate between the chirality 
(1 ± 75) and helicity, which is the projection of the spin of the particle on its momentum. 
In the limit of massless fermions they coincide, but for massive fermions they do not (see 
Problem 1 in the next section). Let us assume that neutrinos are massless and consider 
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the decay 7r+ + ui of pions at rest, where is /i"*" or e"*" and ui is the corresponding 

neutrino. The V — A structure of weak interactions requires the emitted ui to be of left 
handed chirahty. For m,, = this also means that it has the left handed (or negative) helicity 
with its spin antiparallel to its momentum. Conservation of total angular momentum then 
requires Z"*" to be of negative helicity, too (see fig. ||). However, the are antiparticles, 
and due to the V — A structure of weak interactions they must be produced in the states of 
right handed chirality. Therefore the amplitude of the process must be proportional to the 
admixture of right handed chirality in left handed (negative) helicity of the charged lepton, 
i.e. to its mass: A{n'^ — > + ui) oc m;. We therefore expect 

- IK ^"l^ = {-)'{ = 1-28 X 10- . (10) 

where we have taken into account the difference of the phase space factors for the decays 
into and . Experimentally we have 

K = (1.230 ±0.004) X 10-\ (11) 

in a very good agreement with the theoretical prediction; the 4% difference is actually due 
to the fact that (p!0| ) does not include radiation corrections. 

Neutrinos are elusive particles, and neutrino experiments are very difficult. This is 
because the interactions of neutrinos are mediated by heavy and Z° bosons and so at 
low energies they are very weak. The mean free path of a 1 MeV neutrino in lead is about 
1 light year! Therefore neutrino detection requires very large detectors and/or very intense 
neutrino beams. 

What are the main sources of neutrinos? One of the strongest neutrino sources is our 
sun. It emits about 2 x 10^^ electron neutrinos per second, leading to the neutrino flux at 
the surface of the earth of ~ 6 x 10^° cm~'^s~^ in the energy range E < 0.42 MeV and 
~ 5 X 10^ cm~^s~^ in the energy range 0.8 MeV ^ -E < 15 MeV. Nuclear power plants 
are a powerful source of electron antineutrinos. A 3 GW plant emits about 7.7 x 10^° z/g 
of the energy ~ a few MeV per second and creates a flux of ~ 6 x 10^^ cm~^s~^ at 100 
m. At remote locations the averaged number density created by all nuclear power plants 
on the earth is about 10^ — 10'' cm~^. The relic neutrinos, i.e. the neutrinos left over 
from the early epochs of the evolution of the universe, have a number density of about 110 
cm~^ for each neutrino species and a black-body spectrum with the average energy of about 
5 X 10^'* eV. Natural radioactivity of the earth results in the flux of ~ 6 x 10^ cm~^s~^ 
(or number density of ~ 2 x 10^^ cm~^) of neutrinos of the energy E < 1 MeV. The flux 
of atmospheric neutrinos at the earth's surface is ~ 10"^ cm~'^s~^. An important tool for 
studying neutrino properties are particle accelerators which produce neutrino beams of the 
energy ranging typically between 30 MeV and 30 GeV. Finally, a very rare and short-lived 
but yet a very important source of neutrinos are type II supernovae. They emit ~ 6 x 10^^ 
neutrinos and antineutrinos of all flavours over the time interval of about 10 s and with the 
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typical energies E < 30 MeV. Observation of the neutrino burst from the supernova 1987A 
allowed us to obtain important constraints on neutrino properties. 



3 Neutrino mass phenomenology: Weyl, Dirac and 
Majorana neutrinos 

A massless neutrino (or any other fermion) is described by a two-component Weyl spinor 
field. It has a certain chirality, L or R. The chirality projector operators Pl,r are defined 
as 

= P-^-'-^. (12) 

and have the following properties: 

PI = Pl, Pr = Pr, PlPr = PrPl = 0, Pl + Pr = I . (13) 

The terms "left handed" and "right handed" originate from the fact that for relativistic 
particles chirality almost coincides with helicity defined as the projection of the spin of the 
particle on its momentum. The corresponding projection operators are 

They satisfy relations similar to (0). For a free fermion, helicity is conserved but chirality 
in general is not; it is only conserved in the limit m = when it coincides with helicity. 
However, for relativistic particles chirality is nearly conserved and the description in terms 
of chiral states is useful. 

Problem 1. Using plane wave solutions of the Dirac equation show that for positive energy 
solutions, in the limit m = 0, eigenstates of L (R) chirality coincide with the eigenstates of 
negative (positive) helicity. 

For our discussion we will need the particle - antiparticle conjugation operator C. Its 
action on a fermion field ip is defined through 

C: ^^^- = Ci;^ ^ = 1^210. (15) 

The matrix C has the properties 

Ct = = = -C , C^f^C-' = -7j . (16) 

Some useful relations based on these properties are 

{ijy = ij, r = Vc, ^V'2 = ^V'i, j'iAiiJ2 = r2{CA^c~')ri. (ir) 
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where ip, ipi, il)2 are 4-component fermion fields and A is an arbitrary 4x4 matrix. Using 
the commutation properties of the Dirac 7 matrices it is easy to see that, acting on a chiral 
field, C fiips its chirality: 

C: i^L^ {i^iY = mn , i^R ^ H^rT = ir)L , (18) 

i.e. the antiparticle of a left handed fermion is right handed. 

Problem 2. Prove eqs. ^Dj ) and ^T^). 

The particle - antiparticle conjugation operation C must not be confused with the charge 
conjugation operation C which, by definition, fiips all the charge-like quantum numbers of a 
field (electric charge, baryon number, lepton number, etc.) but leaves all the other quantum 
numbers (e.g., chirality) intact. In particular, charge conjugation would take a left handed 
neutrino into a left handed antineutrino that does not exist, which is a consequence of the 
C-noninvariance of weak interactions. At the same time, particle - antiparticle conjugation 
converts a left handed neutrino into right handed antineutrino which does exist and is the 
antiparticle of the left handed neutrino. 

We are now ready to discuss the Dirac and Majorana mass terms. For a massive fermion, 
the mass term in the Lagrangian has the form 

- Cm = mijjip = {ipL + i'R){i'L + i'R) = IpLi^R + '4'R'4'L , (19) 

Thus, the mass terms couple the left handed and right handed components of the fermion 
field, and therefore a massive field must have both components: 

i^ = i^L + i^R. (20) 

Now, there are essentially two possibilities. First, the right handed component of a massive 
field can be completely independent of the left handed one; in this case we have a Dirac 
field. Second, the right handed field can be just a. C - conjugate of the left handed one: 

'4'R = ii'LY = {i''')R, or 

^ = ^L+Vi^lR = ^L+Vi4^LT. (21) 

where we have included the phase factor rj = e**^ with an arbitrary phase (p. In this case 
we have a Majorana field; one can construct it with just one Weyl field. From (|2l|) it 
immediately follows that the C - conjugate field coincides with itself up to a phase factor: 

V'^ = r]*i{j . (22) 

This means that particles described by Majorana fields are genuinely neutral, i.e. coincide 
with their antiparticles. Thus, Majorana particles are fermionic analogs of photons and 7r° 
mesons. To construct a massive Dirac field, one needs two independent 2-component Weyl 
fields, ipL and ipji] together with their C*-conjugates, {ipiY = 4''^r and {ipnY = ip'i, this gives 
four degrees of freedom. In contrast with this, a Majorana fermion has only two degrees of 
freedom, ip^ and {iPlY = V'i?- 
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For massive fermions, particle - antiparticle conjugation C and charge conjugation C 
coincide. For Dirac fermions we have 

c : = ^ ^ = v^L + v^fi = mL + mR, (23) 

where tilde means charge conjugation. For Majorana neutrinos, both particle - antiparticle 
conjugation and charge conjugation leave the field unchanged because it does not have any 
charges 0. As we have already pointed out, particle - antiparticle and charge conjugations 
are not equivalent when acting on chiral fields. 

For n fermion species (fiavours), the Majorana mass term can be written as 



2 



1 



2 



1 
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(24) 

where ip = {ipi, ...,ipn)'^ is a vector in the fiavour space and M is a n x n matrix. Using 
the anticommutation property of the fermion fields and eq. ([T6|), it is easy to show that the 
matrix M must be symmetric: Mij = Mji. It is interesting to note that in classical field 
theory it would be antisymmetric, Mij = —Mji. This, in particular, means that in the case 
of just one neutrino species the Majorana mass vanishes identically in classical field theory. 
It is, therefore, an essentially quantum quantity. Kinematically, Dirac and Majorana masses 
are indistinguishable: they lead to the same relation between energy, momentum and mass 
of the particle, E = v^p^ + m^. 

From eq. (p^) a very important difference between the Dirac and Majorana mass terms 
follows. The Dirac mass terms ipip are invariant with respect to the U{1) transformations 

ip e'^ip , ip ipe-'" , (25) 

i.e. they conserve the corresponding charges (electric charge, lepton or baryon number, etc.). 
It follows from (0) that the Majorana mass terms break all the charges that the field ip has 
by two units. This, in particular, means that, since the electric charge is exactly conserved, 
no charged particle can have Majorana mass. Therefore, out of all known fermions, only 
neutrinos can be Majorana particles. If neutrinos have Majorana masses, the total lepton 
number is not conserved, while it is conserved if neutrinos are Dirac particles. 



4 Neutrino masses in the standard model and slightly 
beyond 

In the standard model of electroweak interactions, all quarks and charged fermions get their 
masses through the Yukawa couplings with the Higgs field H = {H^, H^)^: 

-Cy = hl^QL.UR. H + hi^Q^^R.H + fgL^eR,H + h.c. (26) 

^There may be, however, some differences in the phase factors, see [ p2| . 
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where Qli and Ili are left handed quark and lepton doublets, m^, dRi and cri are SU(2)l - 
singlet right-handed fields of up-type quarks, down-type quarks and charged leptons respec- 
tively, H = iT2H*, Ta being the isospin Pauli matrices, and i, j are the generation indices. 
After the electroweak symmetry is broken by a nonzero vacuum expectation value (VEV) v 
of the Higgs field, the Yukawa terms in ( p6D yield the mass matrices of quarks and charged 
leptons 

irnu)ij = h'ljv , {md)ij = hf-v , {me)ij = f-jV . (27) 

Neutrinos are massless in the minimal standard model. They cannot have Dirac masses 
because there are no SU(2)l - singlet ("sterile") right-handed neutrinos ur in the standard 
model. Can neutrinos have Majorana masses? The answer is no. The reason for this is 
rather subtle, and it is worth discussing it in some detail. 

The Majorana mass term should be of the form uJ^Cul [see eq. (^4])]. Since ui has 
the weak isospin projection I3 = 1/2, the Majorana mass term has I3 = 1, i.e. it is a 
component of the isotriplet operator t^Cir^Tl ~ (3, —2). Therefore, in order to introduce 
the Majorana mass in a gauge invariant way so as not to spoil the renormalizability of the 
standard model, one would need an isotriplet Higgs field A ~ (3, 2): 

- C^^k = fi^ifCiT2Tl)/:^ + h.c. (28) 

Problem 3. Show that for any SU(2) spinorx, 'X^i'T2 has the same transformation properties 
as x\ <md in particular that both X1X2 and x[iT2X2 are invariants while Xi'''X2 and x{iT2TX2 
transform as vectors. Hint: use the property T2T^T2 = — r of the Pauli matrices. 

When the electrically neutral component of A develops a VEV, Majorana neutrino mass is 
generated. However, such a Higgs does not exist in the standard model. Can one construct a 
composite triplet Higgs operator out of two Higgs doublets? Yes, the operator H'^iT2T H ~ 
(3, 2), i.e. has the correct quantum numbers. However, the term {l^C iT2T l){H'^iT2T H) 
has the dimension d = 5, i.e. it cannot enter in the Lagrangian of a renormalizable model 
at the fundamental level. Can the operator 

j^{fClT2Tl){H^lT2TH) (29) 

be generated as an effective operator at some higher loop level? If so, it would produce the 
Majorana mass term for neutrinos ~ fv'^jM when the Higgs field develops a nonva- 
nishing VEV, with M being the characteristic mass scale of the particles in the loop. In 
principle, this is possible. However, in the standard model this does not happen because 
the total lepton number L (more precisely, the difference of the baryon and lepton numbers 
B — L) is exactly conserved. 

Let us discuss this point in more detail. In the standard model, lepton and baryon num- 
bers are conserved at the perturbative level due to accidental symmetries of the Lagrangian. 
These symmetries are called accidental because they are not imposed on the Lagrangian but 
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are just the consequence of the particle content of the standard model, its gauge invariance, 
renormalizability and Lorentz invariance. Indeed, consider, for example, the lepton number 
violating operators II where for simplicity the isospin structure has been suppressed. Such 
operators have isotriplet and isosinglet components: (2, —1) ® (2, —1) ~ (3, —2) © (1, —2), 
both of them with hypercharge —2. Thus, we need a hypercharge 2 field (or product of the 
fields) to obtain a gauge invariant expression. One can, e.g., try ejiU which is the gauge 
invariant AL = 1 operator. However it is a product of three fermion fields and therefore is 
neither Lorentz invariant nor renormalizable. The lepton number violating (AL = 2) oper- 
ator in eq. (^) is both gauge and Lorentz invariant, but, as we have already mentioned, 
it is a c? = 5 nonrenormalizable operator. One can construct, e.g., a baryon number and 
lepton number violating {AB = AL = 1) gauge and Lorentz invariant operator 

£«^^«^C4)(M]^^CeR), (30) 

where a,/3,7 are the SU{3)c indices; however it is a = 6 nonrenormalizable operator. 
Similarly, one can consider all other possible L- and i?-violating operators that can be 
constructed with the standard-model fields and make sure that all of them either violate 
gauge or Lorentz invariance, or are nonrenormalizable. Therefore L and B are automatic 
symmetries of the standard model at the perturbative level. 

Nonperturbatively, both B and L are violated by the electroweak sphalerons since these 
symmetries are anomalous. This violation is tiny and of no practical consequences at low 



energy and temperatures but may be very important, e.g., in the early universe [23|. The 



triangle anomalies of the baryon and lepton number currents are equal to each other in the 
standard model, so that B — L is anomaly free and therefore conserved exactly. Since the 
operators of the type (p9|) which could produce Majorana mass terms for neutrinos break 
not only L but also B — Lhj two units, they cannot be induced in the standard model even 
nonperturbatively. Thus, neutrinos are exactly massless in the minimal standard model. 

It is interesting to note that if the lepton number (or B — L) conservation was not exact, 
neutrino Majorana mass term could be generated even in the standard model framework. 
It has been speculated that quantum gravity should respect only gauge symmetries while 
violating all the global symmetries, such as baryon or lepton numbers. If so, it could induce 
the operators of the type ( |29| ) in the standard model with the mass scale M of the order 
of the Planck mass Mpi ~ 1.2 x 10^^ GeV |24]. This would give neutrino Majorana masses 
m ~ f^/Mpi ~ 10~^ eV, which is exactly of the right order of magnitude to account for the 
solar neutrino deficiency through the vacuum neutrino oscillations. 

How can we extend the standard model so as to accommodate nonzero neutrino mass? 
One can either extend the Higgs content of the model, or the fermion content, or both, or 
enlarge the gauge group (which also requires extended particle content). One possibility 
would be to introduce the isotriplet Higgs field A discussed above. However, such a theory 
would not be very attractive. The VEV u of the triplet Higgs would modify the masses 
of the W and Z bosons at the tree level, and in order not to run into a contradiction 
with the experimental value of the parameter p = / cos^ Oy/m^ it should be < a few 
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GeV. The Yukawa coupling (|28|) conserves the lepton number provided that one assigns the 
lepton charge —2 to the isotriplet Higgs. If the Higgs potential is also L-conserving, the 
VEV u would break the lepton number spontaneously and therefore would give rise to the 



corresponding massless neutral Goldstone boson J which is called "triplet Majoron" ^5 
(in fact, due to the electroweak symmetry breaking, the Majoron is a linear combination 
of the isodublet and isotriplet fields with the dominant triplet contribution). In this case 
even more stringent constraint on the VEV of the triplet Higgs results from astrophysics, 
u < 5 keV. Such small a value of u means that the massive counterpart of the Majoron, p, 
should be very light and can be produced in the decays Z° — > J + p. Such a decay would 
add to the invisible width of the Z° boson a quantity which is equal to the width of decay 
into two extra neutrino species, in sharp contradiction with eq. (0). Therefore the triplet 
Majoron model is ruled out now. If one allows the L-breaking terms in the Higgs potential 
(e.g., p{H^iT2T H)A.^ + /i.e.), lepton number would be broken explicitly and no Majoron 
appears. However, in the standard model, it is highly unnatural to have the triplet Higgs 
VEV u <^v - there is no symmetry or dynamical mechanism to protect it, and its expected 
value is u ~ f , which would lead to a grossly wrong value of the p parameter. 

It is interesting to note that isotriplet Higgs fields may be quite natural in the framework 
of extended gauge theories. For example, the Higgs content of the most popular version of 
the left-right symmetric theories p6| based on the gauge group SU {2)l x SU{2)r x U{1)b-l 
includes a bi-doublet $, which gives the usual Dirac masses mo of fermions, and the fields 
A^, and A/j which are triplets of SU{2)l and SU{2)ji respectively. Their VEVs ul and 
ur produce the Majorana mass terms for the left-handed and right-handed neutrinos, 
and upi. The VEV ur determines the mass scale at which parity is spontaneously broken, 
and is expected to be very large, thereby explaining non-observation of the right-handed 
{V + A) interactions at low energies. Right handed neutrinos are also very heavy in this 



model, while the masses of ul are naturally very small due to the seesaw mechanism [27 



niL m'jj/Mfi (we shall discuss this mechanism in detail in sec. P). The minimization of 
the Higgs potential leads to ul ~ v'^/ur, i.e. the smallness of ul is a natural consequence 
of the largeness of ur, and a kind of the seesaw mechanism for the Higgs VEVs operates. 

Another possible extension of the standard model is just to add three SU{2)l singlet 
neutrinos vr, one per fermion generation. In a sense, this is natural since it restores the 
quark-lepton symmetry: in the minimal standard model all left handed (isodoublet) quarks 
and charged leptons have right handed (isosinglet) counterparts, while neutrinos do not Q 

Right handed neutrinos can have the usual Yukawa couplings to lepton doublets. Since 
they have zero isospin and are electrically neutral, their hypercharge is also zero, i.e. they 
are electroweak singlets. Therefore they can have "bare" Majorana mass terms which are 
invariant with respect to SU{2)i x f/(l)y. The part of the Lagrangian that is relevant for 

■^This argument, however, has to be taken with some caution: the existence of the right handed counter- 
parts of charged left handed fermions is dictated by the requirement of the electroweak anomaly cancellation 
whereas this argument does not apply to right handed neutrinos, see below. 
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the lepton mass generation is therefore 



- £y+M = ftjlueR-jH + iLiPRj H + ^M,, vliCvR, + h.c. (31) 

The second term here yields the Dirac neutrino mass matrix {m£))ij = fJjV. If one assigns 
the lepton charge L = +1 to all the ur fields, the first two terms on the r.h.s. of eq. ( [31| ) 
conserve the total lepton number L [although not the individual lepton numbers, or lepton 
"flavours" Li {i = e,/i, r)]. However, the third term breaks L by two units. Alternatively, 
if one assigns zero lepton charge to i^r, the third term in (|3lD conserves L, but then the 
second term breaks it. Therefore the introduction of ur in the standard model leads to 
a qualitatively new situation: lepton number is no longer an automatic symmetry of the 
Lagrangian following from gauge and Lorentz invariance and renormalizability. 

Since ur are electroweak singlets, they do not contribute to the electroweak anomaly, 
and their number is not fixed by the requirement of the anomaly cancellation. In particular, 
their number need not coincide with the number of generations Ng-. it may as well be smaller 
or even larger than Ng. There are, however, astrophysical and cosmological constraints on 



their number that depend on their masses and mixing parameters ||10| , |Tl|, [12 . 

Problem 4- Assuming that there are two right-handed neutrinos per fermion generation and 
that the Majorana mass terms are absent, write down the most general neutrino mass term 
and analyze it. What would be the neutrino mass spectrum in this case? What would happen 
if there were instead two right-handed charged leptons per fermion generation? 



5 General Dirac + Majorana case 

We shall consider now the most general neutrino mass term for the case of n species of left 
handed and right handed neutrinos. It includes not only the Dirac mass m£, and Majorana 
mass niR for but also the Majorana mass mi for left handed neutrinos. The neutrino 
mass term can be written as 

- Cm = ^1^1 CmLi^L + i^L m}) i^R + ^^^^ C + h.c. = ^^l CM ul + h.c. . (32) 

Here rii = {yi, {urY) = {i^l, is the vector of 2n left handed fields (which we have 
written here as a line rather than column to save space), mi and mR are complex symmetric 
n X n matrices, mo is a complex n x n matrix, and we have introduced the Dirac and right 
handed Majorana mass matrices through their complex conjugates for simplicity of the 
further notation. The matrix has the form 

M = ( 2^ 2- ] , (33) 
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and in deriving eq. (^2]) we have used the relations 



(ul Cm* urY = {v% Cm {v')l , IT^m z/^ = {v')l CmvL = vlCm^ {v'^)^ • (34) 



Problem 5. Prove these relations using properties p^ j and ( |7?[ ) of the matrix C . 

It is instructive to consider first the simple one-generation case in which m/,, mpi and 
m£, are just numbers, and is a 2 x 2 matrix. For simplicity we shall assume all the 
mass parameters to be real. The matrix M. can be diagonalized by the transformation 
M.U = M-d where U is an orthogonal matrix and J^d = diag{mi, m2). We introduce 
the fields Xl through rii = Uxl, or 

J \ -smO COSO J \ X2L J 

Here Xil and X2L are the left handed components of neutrino mass eigenstates. The mixing 
angle 6 is given by 

tan 2^ = 36 

mR-mi 

and the neutrino mass eigenvalues are 



mi^2 = T \ +mj^. (37) 



They are real but can be of either sign. The mass term can now be rewritten as 

- Cm = ^nlCMnL + h.c. = ^XLCMdXL + h.c. 

1 1 
= 2("^i^i'^'^^i^ + "^2xIlC'X2l) + /i-c. = -( |mi| XiXi + l"^2|X2X2) • (38) 

Here we have defined 



Xl = XiL + ViiXiiY , X2 = X2L + V2{X2Ly ■ (39) 

with ?7i 2 = 1 or — 1 for mi 2 > or < respectively. It follows immediately from eq. ( PPj ) 
that the mass eigenstates xi and X2 are Majorana neutrinos. The relative signs of the mass 
eigenvalues (771 and 772) determine the relative CP parities of xi and X2', physical masses |mi| 
and |m2| are positive, as they should. The fact that the mass eigenstates in the case of the 
most general Dirac + Majorana mass term are Majorana particles should not be surprising. 
We have four neutrino degrees of freedom, u^, and their C- conjugates {i^lY = ^'r 
and {i^rY = v'^i- The mass matrix has two different eigenvalues, so there are two massive 
neutrino fields. Each of them corresponds to two degrees of freedom, therefore they have 
to be Majorana particles. Analogously, in the case of n generations, the most general mass 
term (|3^ ) leads to 2n massive Majorana neutrinos. 
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It is instructive to consider some limiting cases. In the limit of no Majorana masses 
{rriL = fTLji = 0), pure Dirac case has to be recovered. In this case the mass matrix ( P^D has 
the form 

This matrix corresponds to a conserved lepton number Lj^^ — L^c^ = L^^ + Lj,^ which can 
be identified with the total lepton number L. Thus, the lepton number is conserved in this 
limitine; expected. Let us now check that the usual Dirac mass term is recovered. 

The matrix M. in P0| ) is diagonalized by the rotation (|35| ) with 9 = 45°, and its eigen- 
values are —m and m. This means that we have two Majorana neutrinos that have the same 
mass, opposite CP parities and are maximally mixed. We now want to demonstrate that this 
is equivalent to having one Dirac neutrino of mass m. We have ri2 = —rji = 1; from eqs. (|35| ) 
and (^) it then follows xi + X2 = V2(z/l + ur), xi - X2 = -V2{iyl + z/'^j) = -(^i + 
This gives 

^ (XiXi + X2X2) = ^"^[(Xi + X2)(Xi + X2) + [(Xi -X2)(xi -X2)] =muDi^D, (41) 
where 

ud = ul + i^r- (42) 

The counting of the degrees of freedom also shows that we must have a Dirac neutrino 
in this case - there are four degrees of freedom and just one physical mass. Thus, two 
maximally mixed degenerate Majorana neutrinos of opposite CP parities merge to form a 
Dirac neutrino. It is easy to see that in this case their contributions into the probability 
amplitude of the neutrinoless double beta decay 2/30z/ (which is only possible for Majorana 
neutrinos) exactly cancel, see eq. ( |126|) in sec. ||. 



If the Majorana mass parameters rriL and m^j do not vanish but are small compared to 
m£), the resulting pair of Majorana neutrinos will be quasi-degenerate with almost maximal 
mixing and opposite CP parities. The physical neutrino masses in this case are |?7ii 2| — 
m^) ± {rriL + m^)/2. Such a pair in many respects behaves as a Dirac neutrino and therefore 
sometimes is called a "quasi-Dirac neutrino". In particular, its contribution to the 2/50z/ 
decay amplitude (|126|) is proportional to the mass difference (m^ + mR)/2 which is much 



smaller than the mass of each component. 



6 Seesaw mechanism of neutrino mass generation 

The natural mass scale in the standard model is the electroweak scale which is of the order 
of f ^ 174 GeV. The smallness of, e.g., the electron mass rrie — 0.511 MeV is not explained; 
however, it is easily accommodated in the standard model through the proper choice of the 
corresponding Yukawa coupling, /g ~ 3 x 10~^. At the same time, similar explanation of 
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the smallness of the electron neutrino mass, m,y^ < 5 eV, would require the Yukawa coupling 
10~^^. Does this pose any problem? If we are willing to accept a very small Yukawa 
coupling of the electron, why should not we accept small neutrino Yukawa couplings as well? 
After all, 10~^^ may be as good (or as bad) as 10~^. 

The problem is that, except for neutrinos, the masses of all the fermions in each of 
the three generations are within 1-2 orders of magnitude of each other. The inclusion of 
neutrinos leads to huge disparities of the fermion masses within each generation. Therefore, 
if a future more complete theory explains why there is a large mass hierarchy between 
generations, it would still remain to be explained why neutrinos are so light compared to 
the other fermions of the same generation. 

The seesaw mechanism |^ provides a very simple and attractive explanation of the 
smallness of neutrino mass. It relates the smallness of my with the existence of very large 
mass scales. Although the seesaw mechanism is most natural in the framework of the grand 
unified theories (such as 50 (10)) or left-right symmetric models, it also operates in the 
standard model extended to include the right handed ("sterile") neutrinos vr. This is the 
model discussed in the end of sec. |[ and the relevant part of the Lagrangian is given by 
eq. (pri). We shall now consider the seesaw mechanism in some detail. 

The most general mass term for n generations of left handed and right handed neutrinos 
is written in eq. (|32[) with the mass matrix given in eq. (|33|) . Notice that in the standard 
model there is no Majorana mass term for left handed neutrinos since there are no triplet 
Higgs scalars; however, is different from zero in some extensions of the standard model, so 
we shall keep it for generality. The right handed neutrino vr is an electroweak singlet and so 
its mass is not protected by the electroweak symmetry. One can therefore expect it to be very 
large, possibly at the Planck scale or at some intermediate scale Mj ~ ^JvM■p\ ~ 10^° — 10^^ 
GeV which may be relevant for the physics of parity breaking. 

Let us first consider the limit m^ <^ <^ m^ of the simple one-generation case 
discussed in the previous section. In this limit 

9 ^ ^ 1 , mi ~ mi — , m2 — m^ , (43) 

rriR mR 

xi-^L + mi^iT , X2 ^ {^rY + m^R ■ (44) 

Thus we have a very light Majorana mass eigenstate xi predominantly composed of ul and 
a heavy eigenstate X2 mainly composed of ur. The admixture of the singlet neutrino state 
ur in xi and that of the usual neutrinos ul in X2 are of the order of mD/mR. As follows 
from eq. (^Sp, it is ur being heavy that makes light 

■^This is the commonly used jargon: in fact, i/^ and vr, being chiral fields, do not have any mass but 
rather are components of some massive fields, ft would be more correct to say that the fact that X2 made 
predominantly of vr and its C-conjugate {vf!.Y heavy explains the lightness of xi made predominantly 
of vl and {i^lY, but this is just too lengthy. 
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Consider now the full n-generation case. We want to block-diagonalize the matrix Ai 
in eq. (^Sp so as to decouple light and heavy neutrino degrees of freedom: 

where ?7 is a unitary 2n x 2n matrix, and we have changed the notation — >■ Mn- We 
shall be looking for the matrix U of the form 

f/= ( J ) , f/tf/ = 1 + , (46) 

where the elements are n x n matrices, and p will be treated as a perturbation. We shall 
neglect for simplicity possible CP violation in the leptonic sector and take m^, rriD and Mji 
to be real matrices (effects of CP violation in neutrino oscillations will be discussed in sec. 
|7.3|) . The matrix p can then also be chosen to be real. Block-diagonalization of Ai gives 

p ~ tudM^^ , friL^niL- mDM]^^m^ , MrC:^ Mr. (47) 

These relations generalize those of eq. (^31) to the case of n generations. The diagonalization 
of the effective mass matrix fhi yields n light Majorana neutrinos which are predominantly 
composed of the usual ("active") neutrinos vl with very small (~ mD/MR) admixture of 
"sterile" neutrinos ur; diagonalization of Mr produces n heavy Majorana neutrinos which 
are mainly composed of ur. It is important that the active neutrinos get Majorana masses 
rfiL even if they have no "direct" masses, i.e. rriL = 0, as it is in the standard model. The 
masses of active neutrinos are then of the order of mj^/MR. Generation of the effective 
Majorana mass of light neutrinos is diagrammatically illustrated in fig. |^. It is interesting 
that with the largest Dirac mass eigenvalue of the order of the electroweak scale, m^j ~ 200 
GeV, the right handed scale Mr ~ 10^^ GeV which is close to the typical GUT scales, and 
assuming that the direct mass term < mj^/MR, one obtains the mass of the heaviest of 
the light neutrinos ~ (10~^ — 10~^) eV, which is just of the right order of magnitude for 
the neutrino oscillation solution of the atmospheric neutrino anomaly. 

(H) {H) 



» 1 » X < 1 < 

UL rUD l^R Mr mo I'L 

Figure 3: Seesaw mechanism of niL generation 

Problem 6. Perform the approximate block diagonalization of the matrix Ai and verify 
eq. 
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As was pointed out before, in general the number k of "sterile" neutrino species need 
not coincide with the number n of the "active" ones. The seesaw mechanism works when 
7^ 77, as well, and the formulas of eqs. (|45|) - (|47| ) are still valid. The difference is that 
now the matrices m^) and p are rectangular nx k matrices (n lines, k columns) rather than 
square matrices. The matrices mi and Mr are square matrices of dimension nxn and kx k 
respectively; the same is true for the unit matrices in the left upper and right lower corners 
of the matrix U in (^61). 

In deriving eq. we have assumed that the matrix p can be considered as a small 
parameter, i.e. that Mji ^ mjj. In what sense one can consider one matrix to be much 
larger than another? Does that mean that we have to require that all the elements of 
be much larger than all the elements of m/)? Obviously, this is not correct. One can, for 
example, consider the matrix Ai in the basis of in which Mr has been diagonalized; in 
this case all off-diagonal elements of Mji are zero and yet the seesaw approximation works 
perfectly well. On the contrary, if one chooses the matrix Mji with all its elements equal to 
each other, the standard seesaw mechanism fails, no matter how large the elements of Mr. 
This is just because det{MR) = in this case, so that the matrix M^^ does not exist. The 
true criterion of the applicability of the seesaw approximation is that all eigenvalues of Mr 
be large compared to all eigenvalues of m£,. 

Problem 7. Consider the seesaw mechanism for three ul and three vr species in the case 
when one of the eigenvalues of Mr is zero. What will he the neutrino mass spectrum in this 
case? Hint: go to the vr basis where Mr is diagonal and include the line and the column 
that contain zero eigenvalue into a redefined matrix thl. 

7 Neutrino oscillations in vacuum 

Neutrino oscillations are the most sensitive probe of neutrino mass. Solar and supernova 
neutrino experiments may be able to discover the neutrino mass as small as 10~^ eV or 
even smaller, far beyond the reach of the direct kinematic search experiments. The idea 
of neutrino oscillations was first introduced by Pontecorvo [^. The essence of this effect 
is very simple, and examples of the oscillation phenomena can be found in nearly every 
textbook on quantum mechanics. Consider, for example, a two-level quantum system. If 
the system is in one of its stationary states (eigenstates of the Hamiltonian) , it will 
remain in this state, and the time evolution of the wave function is very simple - it just 
picks up a phase factor: \^!i{t)) = e~*'^'* 1^^,(0)). If, however, a state is prepared which is not 
one of the eigenstates of the Hamiltonian of the system, the probability to find the system 
in this state will oscillate in time with the frequency = E2 — Ei where Ei and E2 are 
the eigenenergies of the system. 

In the case of neutrino oscillations, neutrinos are produced by the charged-current weak 
interactions and therefore are weak-eigenstate neutrinos Ve-, or v^-. However, the neutrino 
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mass matrix in this (flavour) basis is in general not diagonal. This means that the mass 
eigenstate neutrinos i>i, z/2 and z/3 (the states that diagonalize the neutrino mass matrix, 
i.e. the free propagation eigenstates) are in general different from the flavour eigenstates. 
Therefore the probability of finding a neutrino created in a given flavour state to be in the 
same state (or any other flavour state) oscillates with time. 

We shall consider neutrino oscillations in the case of Dirac neutrino mass term and then 
will comment on the Majorana and Dirac + Majorana cases. The part of the Lagrangian 
that describes the lepton masses and charged current interactions is 

- Cw+m = ^ <L Y ^'aL + i^l) abe'aLe'bR + {mD)ab KlKr + ^.C. (48) 

Here the primes are used to denote the flavour eigenstate fields. It follows from this ex- 
pression that the individual lepton flavours Le, and are not conserved when the 
Dirac neutrino mass term is present, but the total lepton number L = Lg + + L^- is still 
conserved. 

The mass matrix of charged leptons mi and the neutrino mass matrix itid are general 
complex matrices which can be diagonalized by bi-unitary transformations. Let us write 



e 



'l = VLeL, = Vr cr v'^ = UlI'l, i^'r = Ur vr , (49) 



and choose the unitary matrices Vl, Vr, Ul and Ur so that they diagonalize the mass 
matrices of charged leptons and neutrinos: 

Vl nil Vr = {mi)diag , ul rUD Ur = {mD)diag ■ (50) 

The "unprimed" fields en, CiR, and uir are then the components of the Dirac mass 
eigenstate fields e, = en + eiR and Vi = i^iL + ^iR- The Lagrangian in eq. (^8[) can be written 
in the mass eigenstate basis as 

- Cw+rn = (ylUL)ij I'Lj W~ + mueLieRi + rriDi VnVRi + h.c. , (51) 

where mu are the charged lepton masses and m£,i are the neutrino masses. The matrix 
U = V^Ul is called the lepton mixing matrix, or Maki-Nakagawa-Sakata (MNS) matrix 



29| . It is the leptonic analog of the CKM mixing matrix. It relates a neutrino flavour 
eigenstate produced or absorbed alongside with the corresponding charged lepton, to 
the mass eigenstates 

W'a) = U:,W^), (52) 

In what follows, to simplify the notation we shall omit the primes and distinguish between 
the flavour and mass eigenstates just by using the indices from the beginning of the Latin 
alphabet for the former and from its middle for the latter. Assume that at a time t = 
the flavour eigenstate |z/a) was produced. What is the probability to find the neutrino in a 
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state \vb) at a later time tl Time evolution of the flavour eigenstates is not simple, therefore 
it is more convenient to follow the evolution of the system in the mass eigenstate basis 
(we shall discuss evolution in the flavour basis in sec. H). The initial state at t = is 
|z/(0)) = |z/a) = Ulj\vj)] the neutrino state at a later time t is then 

kW)=f/;-e-^^^V,). (53) 

The probability amplitude of flnding the neutrino at the time t in a flavour state \vh) is 

A{ua t) = {uMt)) = e-*^^* {ub\u,) = UuKj e-'"^' = Ut, e"^^^* f/^- . (54) 

As usual, the sum over all intermediate states j is implied. The last expression here has a 
very simple physical meaning. The factor U*j = Uj^ is the amplitude of transformation of 
the initial flavour eigenstate neutrino i^a into a mass eigenstate one Uj] the factor e"*'^-'* is 
just the propagator describing the time evolution of the mass eigenstate neutrino Uj in the 
energy representation, and flnally the factor Ubj converts the time-evolved mass eigenstate 
into the flavour eigenstate Ub. The neutrino oscillation probability, i.e. the probability of 
the transformation of a flavour eigenstate neutrino z/^ into another one Ub, is then 

P{u, ^ Ub, t) = \A{u, Ub, t)\' = |f/,,e-*^^*f/:/. (55) 
We have discussed neutrino oscillations in the case of Dirac neutrinos. What happens 



if neutrinos have a Majorana mass term rather than the Dirac one? Eq. (|48|) now has to 
be modifled: the term {jnD)abT^'aL^bR + ^-C- has to be replaced by {'mM)abi^aLKR + — 
{nT'M)ab^Z.C z/^^ + h.c. This mass term breaks not only the individual lepton flavours but 
also the total lepton number. The symmetric Majorana mass matrix {mM)ab is diagonalized 
by the transformation UJ^uim Ul = {rnM)diag, so one can again use the fleld transformations 
(|49|). Therefore the structure of the charged current interactions is the same as in the case of 
the Dirac neutrinos, and the diagonalization of the neutrino mass matrix in the case of the n 
fermion generations again gives n mass eigenstates. Thus the oscillation probabilities in the 
case of the Majorana mass term are the same as in the case of the Dirac mass term. This, 
in particular, means that one cannot distinguish between Dirac and Majorana neutrinos by 
studying neutrino oscillations. Essentially, this is because the total lepton number is not 
violated by the neutrino flavour oscillations. 

The situation is different in the case when the neutrino mass term is of the most general 
Dirac + Majorana (D + M) form (which, like in the pure Dirac case, requires the existence of 
the electroweak singlet neutrinos upt). In particular, in the case of n and n upt species, the 
neutrino mass matrix has the dimension 2n x 2n, leading to 2n massive Majorana neutrino 
states. The total lepton number conservation is violated by the Majorana mass term. Unlike 
in the pure Dirac and pure Majorana cases, in this D + M case one can have a new type of 
neutrino oscillations: in addition to the usual flavour oscillations z/^ —>■ Ub, oscillations into 
"sterile" states, — »■ ^^6, can occur. These oscillations violate the total lepton number L. 
Thus, the D + M case can in principle be distinguished from the pure Dirac and Majorana 
cases in neutrino oscillations experiments. For more detailed discussion of the D + M case 
see, e.g., p 
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7.1 2 flavour case 



Let us now consider neutrino oscillations in a simple case of just two neutrino species, 
and z/^. The lepton mixing matrix U can be written as 

where c = cos^o? s = sin^o? being the mixing angle. The neutrino mass and flavour 
eigenstates are therefore related through 

We) = C\ui)+S\u2), 

=-s\ui) + C\U2) ■ (57) 

Substituting ( pBD into ( ^3] ) and taking into account that for relativistic neutrinos of the 
momentum p, 

2 2 

mf mf 



Ei = ^p"^ + mf c^p+ ^:^p+ (58) 
we find the transition probabilities 

CAm^ \ 
-J^A ■ (59) 

Here Am^ = m\ — m\. The survival probabilities are P{ve —>■ Ve] t) = P{i'^ —>■ i^n', t) = 
1 — P(z/e —>■ Ve'i t). It is convenient to rewrite the transition probability in terms of the 
distance L travelled by neutrinos. For relativistic neutrinos L c::^ t, and one has 

,2 ofl „;„2 ^ L 



P{ue ^ i^/.; L) = sin' 2^0 sin' ( tt— ) , (60) 



''OSC 



where lose is the oscillation length defined as 

A-kE ^ ^„ E(MeV) E(GeV) 

lose = -. — - ~ 2.48 m — ^ = 2.48 km ^ \ ^ — ^ (61) 

Am2 Am2(eV') Am2 (eV') ^ ^ 

It is equal to the distance between any two closest minima or maxima of the transition 
probability (see fig. Notice that lose is inversely proportional to the energy difference 
of the neutrino mass eigenstates: lose = 27r/(i^^2 ~ -^'i)- Another convenient form of the 
expression for the transition probability is 

P(z/e ^ z/^; L) = sin' 200 sin' (^1.27Am'^j , (62) 
where L is in m and E in MeV or L is in km and E in GeV. 
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Figure 4: Pattern of two-flavor neutrino oscillations (neutrino energy uj, distance z). From [13| 



Let us discuss the probability of neutrino oscillations ( pOD (see fig. |[). It has two 
factors. The first one (sin^ 26*0) does not depend on the distance travelled by neutrinos; it 
describes the amplitude (or depth) of the neutrino oscillations. The amplitude is maximal 
(equal to one) when the mixing angle Oq = 45°, which corresponds to maximal mixing. 
When 9o is close to zero or 90°, flavour eigenstates are nearly aligned with mass eigenstates, 
which corresponds to small mixing. In this case the oscillation amplitude is small. The 
second factor oscillates with time or distance L travelled by neutrinos. The phase of the 
sine (the oscillation phase) is proportional to the energy difference of the mass eigenstates 
Am? /2E and to the distance L. In order to have an appreciable transition probability, 
it is not enough to have large mixing: in addition, the oscillation phase should be not 
too small. When the oscillation phase is very large, the transition probability undergoes 
fast oscillations. Averaging over small energy intervals corresponding to the finite energy 
resolution of the detector, or over small variations of the distance between the neutrino 
production and detection points corresponding to the finite sizes of the neutrino source and 
detector, results then in averaging out the neutrino oscillations. The observed transition 
probability in this case is 



i sin^ Id, 
2 



• 



(63) 



In our discussion we have been assuming that the initially produced neutrino state 
has a certain momentum. Since it is a mixture of different mass eigenstates, its energy is 
not well defined in this case. We then consider the time evolution of the neutrino state. 
Alternatively, one could assume that the neutrino is produced in a state of a certain energy 
but not well defined momentum, and consider the evolution of the neutrino system in space. 
For relativistic neutrinos the result will be the same. Which description is more correct? 
Does the initially produced state have a certain energy or a certain momentum? In general, 
since the production process takes finite time and is localized in space, neither neutrino 
energy nor momentum is well defined; moreover, the transition probability depends on both 
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the neutrino production and detection processes. Fortunately, for relativistic neutrinos these 
subtleties are unimportant, and the simple description in terms of the states with certain 
momentum or energy gives the correct answer 0. 



7.2 3 flavour case 



Consider now the case of three neutrino flavours. The neutrino flavour eigenstate and mass 
eigenstate fields are related through 



(64) 



In general, in the case of Dirac neutrinos the lepton mixing matrix U in (0) depends on 
three mixing angles 612, ^13 and 623 and one CP- violating phase 6 (in the case of Majorana 
neutrinos there are two additional, so-called Majorana phases, see discussions in sections 
773| and p.4|). It is convenient to use the parametrization of the matrix U which coincides 
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with the standard parametrization of the quark mixing matrix ||2T 
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(65) 



Here Cjj = cos 6ij, Sij = sin 6ij. The probabilities of oscillations between various flavour 
states are given by the general expression (0). Unlike in the two-flavour case, they in general 
do not have a simple form. There are, however, several practically important limiting cases in 
which one can obtain very simple approximate expressions for the oscillation probabilities 
in terms of the 2-flavour ones. Assume first that the neutrino mass squared differences 
Ami 



m~ — rm have a hierarchy 



|Am2i| 



|Am32| 



(66) 



This means that either mi <^ (<) m2 <^ ^3 (direct hierarchy) or 7713 <C mi ^ m2 (inverted 
mass hierarchy). These cases are of practical interest since the solar neutrino data indicate 
that one needs a small mass squared difference Attt-q ~ 10~^ eV^ for the solution of the solar 
neutrino problem through the matter-enhanced neutrino oscillations (or Am^ ~ 10"^'^ eV^ 
for the vacuum oscillations solution) whereas the explanation of the atmospheric neutrino 



experiments through the neutrino oscillations requires Am' 



atm 



10 eV , much larger than 



Am?,. Consider first the oscillations over the baselines L for which 



"0- 



^"^21 r 

—L < 1 

2E 



(67) 



^Uncertainties of neutrino energy and momentum may be important for the question of coherence of 
the neutrino state in the case of the oscillations over very long baselines L, e.g. for oscillations of solar or 
supernova neutrinos on their way to the earth. For discussion see, e.g., [pO[. 



23 



This case is relevant for atmospheric, reactor and accelerator neutrino experiments. It 
follows from ( pTj) that the oscillations due to the small mass difference Am^i are effectively 
frozen in this case, and one can consider the limit Am^i —>■ 0. The probability of z/q —>■ Ub 
oscillations then takes a very simple form: 

P{i^a Vb\ L) = 4|f/,3r lUbsl'' sin^ I ^#^1 ■ (68) 



Ue, Uf, and 
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with Piyh 





h,? \U,a? sin^ = 4 sin^ 2^13 sin^ ^) , (69) 

W \Ur3\' sin^ = 4 sin^ 29,, sin^ , (70) 

If/.ar sin^ ^) = 4 sin^ 29,, sin^ l) , (71) 



of the lepton mixing matrix U and one mass squared difference. The survival probability 
for electron neutrinos takes a particularly simple form 

P{u, ^ z/e; L) = 1 - sin^ 29,, sin' , (72) 

i.e. it coincides with the i^e survival probability in the 2-flavour case with the mass squared 
difference Am^ = Amg^^ and mixing angle 9q = 9,,. 

Consider now another limiting case, which is relevant for the solar neutrino oscillations 
and also for very long baseline reactor experiments (such as KamLAND, see sec. ^2]). We 
shall be again assuming the hierarchy (|66D and in addition 



Amoi AmL , , 

whereas the condition (|67D is no longer necessary. In this case the oscillations due to the 
mass squared differences Am,, and Am32 are very fast and lead to an averaged effect; the 
z/g survival probability is then 

P{u, -> u,) ^ 4,P + 4 , (74) 



where P is the i^e survival probability in the 2-fiavour case with the mass squared difference 

^21 



Am^ = Am2i and mixing angle 9q = 9,2. In the case of neutrino oscillations in vacuum one 



has 



(Ami 



P = 1 - sin^ 2^12 sin^ — L . (75) 
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Finally, consider the limit Ues = (the results will be also approximately valid for If/gsl <^ !)• 
In this case one obtains 

P(i/e — > Uf^; L) = C23 sin^ 26i2 sin^ A21 , (76) 
P(z/e Ur] L) = S23 sin^ 26112 sin^ A21 , (77) 
P(z/^ Ur] L) = sin^ 26123 {s^ c?2 sin^ A21 + sf^ sin^ A31 + cf^ sin^ A32) , (78) 

where Ajj = {Amfj/4E)L, and no assumption about the hierarchy of the mass squared 
differences has been made. Notice that the limiting cases that we have discussed here are 
not mutually excluding, i.e. have some overlap with each other. 

In general, when considering the propagation of solar neutrinos in the sun or in the 
earth, one should take into account matter effects on neutrino oscillations. The same is true 
for the terrestrial atmospheric and long baseline accelerator neutrino oscillation experiments 
in which the neutrino trajectories or their significant portions go through the matter of the 
earth. Matter effects on z/^ ^ Ur oscillations are relatively small (they vanish in the 2-flavour 
approximation), but they may be quite appreciable for z/g ^ and z/g ^ Ur oscillations. 
We shall discuss these effects in sec. |^. As we shall see, eq. (^) remains valid in the case 
of neutrino oscillations in matter as well, but the two-flavour probability ( ffSD has to be 
modified. Eqs. ( |69D - ([7TD, (p2[) and (^) - ( [78| ) are also modified in the case of neutrino 
oscillations in matter (see sec. 



Problem 8. Assuming that conditions ^ddj) and ^dlj) are satisfied derive eq. ^6^). 
Problem 9. Derive eqs. ^7^ ) - ^Tdj ) assuming ^6^ ) and j^T^). 



Problem 10. Assuming U^z = derive eqs. ( \7di ) - (\7Bi). 

Problem 11. Imagine a world described by the standard model supplemented by three right 
handed neutrinos vr, but in which the charged leptons are massless (Yukawa couplings 



in eq. vanish). Will neutrinos oscillate in such a world? 

Problem 12. Do neutrinos produced in the decay — > vv oscillate? (When you come to 
an answer, compare it with that in 



7.3 Implications of CP, T and CPT symmetries 

We shall now consider the consequences of CP, T and CPT symmetries for neutrino os- 
cillations. As we discussed in sec. H, charge conjugation operation C is not well defined 
for neutrinos as it would convert a left handed neutrino into a non-existent left handed 
antineutrino. CP, however, is well defined: it converts a left handed neutrino into a right 
handed antineutrino which is the antiparticle of z/^. Thus, CP essentially acts as the par- 
ticle - antiparticle conjugation. If CP is conserved, the probabilities of oscillations between 
particles and their antiparticles coincide: 

CP : P{Va ^ t) = P{Va ^ i^b\ t) . (79) 



25 



In quantum field theory, the same field operator that annihilates a particle also creates its 
antiparticle, whereas its Hermitean conjugate does the opposite. Therefore the action of 
the particle - antiparticle conjugation on the lepton mixing matrix U amounts to U —>■ U*. 
This means that CP is only conserved in the leptonic sector if the mixing matrix U is real 
or can be made real by a rephasing of the lepton fields. 

In general, a unitary n x n matrix depends on n{n — l)/2 angles and n{n + l)/2 phases. 
In the Dirac case, 2n — 1 phases can be removed by a proper rephasing of the left handed 
fields, leaving n{n + l)/2 — {2n — 1) = {n — l){n — 2)/2 physical phases (the rephasing 
of the left handed field leaves the lepton mass terms unchanged since the phases can be 
absorbed into the corresponding rephasing of the right-handed fields). Thus, in the Dirac 
case CP non-conservation is only possible in the case of n > 3 generations. In the Majorana 
case there is less freedom to rephase the fields since the Majorana mass terms are of the 
form u^Ul + h.c. rather than of the form i'rI'l + h.c. and so the phases of neutrino fields 
cannot be absorbed. Therefore in the Majorana case only n phases can be removed, leaving 
n(n + l)/2 — n = n(n — l)/2 physical phases. Out of these phases, (n — l)(n — 2)/2 are the 
usual, Dirac-type phases while the remaining n — 1 are specific for the Majorana case, so 
called Majorana phases. Majorana phases do not lead to any observable effects for neutrino 
oscillations [1321, and we shall not consider them here. 



CPT transformation can be considered as a combined action of CP, which results in 
the interchange between paricles and antiparticles, and time reversal, which interchanges 
the initial and final states. Therefore under CPT transformation the oscillation probability 
P(fa — Ub] t) goes into P(z/b — > Ua] t). One the other hand, from eq. (p^ it follows that the 
CPT transformation U U* , t —t transforms the oscillation amplitude A(z/a Ub, t) 
into its complex conjugate. Therefore the oscillation probabilities are invariant with respect 
to CPT, i.e. the following equality holds: 

CPT : P{ua ^ i^b] t) = P{ub ^ ya\ t) . (80) 

CPT invariance implies, in particular, that the survival probabilities for neutrinos and an- 
tineutrinos are the same: Piya — ^ T-'a] t) = P{pa ^a', t). 

Finally, time reversal T interchanges the initial and final states, so if T is conserved one 

has 

T : P(z/, ^ z/fe; t) = P{vb ^ Va\ t) . (81) 

From CPT invariance it follows that CP conservation is equivalent to T conservation; indeed, 
it is easy to see that ([79|) and ( pO]) lead to (|81|). 



Let us now consider effects of CP-violation on neutrino oscillations. If CP is not con- 
served, the oscillation probabilities for neutrinos are different from those for antineutrinos. 
As we pointed out before, this is only possible if the lepton mixing matrix U is essentially 
complex, i.e. has unremovable phases. For three lepton generations, there is only one such 
phase S, and so there should be only one CP-odd oscillation asymmetry. Let us denote the 
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CP-odd asymmetries as 

APab = P{l^a ^ l^b, t) - P{Ua ^ ^b] t) . (82) 



From CPT invariance one finds AP^b = —^Pta- Using the parametrization ( p5| ) of tlie 
mixing matrix U it is not difficult to find 

APe^ = AP^r = APre = 4si2 Cu S13 c^g S23 C23 sin 6 



/Am?, \ f AmL \ ( Aml^ ^ 

sin — -^t + sm — -^t + sm — #t 



3) 



2E ) \ 2E ) \ 2E 

Problem 13. Derive eq. using parametrization (1^^ and unitarity property of the lepton 



mixing matrix U . 

This expression has several interesting features. First, as expected, it vanishes in the limit 
(5 = 0. Second, it vanishes if any of the mixing angles ^12, ^13 or ^^23 is zero or 90°. In 
particular, oscillation probabilities (|76D-([78D obtained in the limit ^13 = are CP-symmetric. 
Third, since the mass squared differences satisfy the relation Amfj + Am^g + Am\^ = 0, the 
CP-odd asymmetry ( pSf ) vanishes if even one of Amfj is zero. We have already encountered 
this situation - the transition probabilities (^) - ([72|). derived in the limit Am^i = depend 
only on the absolute values of Uai and are therefore CP-invariant. The relation between Amfj 
also means that in the limit of the small arguments of the sines (small oscillation phases) in 
eq. ( |83D the probability asymmetry is cubic (more precisely, tri-linear) in small phases rather 
than linear. Another important difference between the usual oscillation probabilities and 
the CP-odd asymmetry (P^D is that while the former contain squared sines of the oscillation 
phases and therefore oscillate near non-zero average values, the latter is linear in these sines 
and so oscillates around zero. This means that in the case of very large oscillation phases, 
when the averaging regime sets in, the CP-odd asymmetry of neutrino oscillations averages 
to zero. It is for this reason that the oscillation probability ([7^) is CP-invariant. 



It follows from the above discussion that the experimental observation of CP violation 
effect in neutrino oscillations is a very difficult task. The CP-odd probability asymmetry 
is suppressed if any one of the three lepton mixing angles is small, and we know from the 
CHOOZ reactor neutrino experiment that the mixing angle 613 is small (see sec. ^ below 
for a discussion of the reactor and accelerator data). The atmospheric neutrino experiments 
indicate that the mixing angle ^23 is rather large, but the value of 612 is largely unknown: 
the present solar neutrino data allow both large and small values for this mixing angle. The 
hierarchy ( |66D of the mass squared differences which follows from the solar and atmospheric 
neutrino observations further hinders experimental searches of CP violation in the neutrino 
oscillation experiments. In addition, matter effects on neutrino oscillations may mimic CP 
violation and so make the searches of the genuine CP violation even more difficult (see sec. 
Still, the presently allowed values of neutrino parameters do not exclude a possibility 



of observation of CP violation effects in the future neutrino experiments. Discovering CP- 
nonconservation in the leptonic sector would be of great importance, so the goal is worth 
pursuing. 
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Problem 14- Using CPT invariance and the unitarity condition P{ve ^e) + -P(^'e 
i/fj) + P{Vf, i>T-) = 1 (and similar conditions for and v^-) show that there is only one 
independent CP -odd oscillation asymmetry in the case of three neutrino flavours. 

Problem 15. Apply the same approach to the four-flavour case and find the number of 
independent CP-odd oscillation asymmetries in that case. Compare it with the number of 
the physical Dirac-type phases in the mixing matrix U in the case of four generations. 



8 Neutrino oscillations in matter 

Neutrino oscillations in matter may differ from the oscillations in vacuum in a very significant 
way. The most striking manifestation of the matter effects on neutrino oscillations is the 
resonance enhancement of the oscillation probability - the Mikheyev - Smirnov - Wolfenstein 
(MSW) effect [|, In vacuum, the oscillation probability cannot exceed sin^26'o, and 



for small mixing angles it is always small. Matter can enhance neutrino mixing, and the 
probabilities of neutrino oscillations in matter can be large (close to unity) even if the mixing 
angle in vacuum is very small. Matter enhanced neutrino oscillations provide a very elegant 
solution of the solar neutrino problem; matter effects on oscillations of solar and atmospheric 
neutrinos inside the earth can be quite important. Neutrino oscillations in supernovae and 
in the early universe may also be strongly affected by matter. 

How does the matter affect neutrino propagation? Neutrinos can be absorbed by the 
matter constituents, or scattered off them, changing their momentum and energy. However 
the probabilities of these processes, being proportional to the square of the Fermi constant 
Gp, are typically very small. Neutrinos can also experience forward scattering, an elastic 
scattering in which their momentum is not changed. This process is coherent, and it creates 
mean potentials Va for neutrinos which are proportional to the number densities of the 
scatterers. These potentials are of the first order in Gp, but one could still expect them to be 
too small and of no practical interest. This expectation, however, would be wrong. To assess 
the importance of matter effects on neutrino oscillations, one has to compare the matter- 
induced potentials of neutrinos Va with the characteristic neutrino kinetic energy differences 
Am'^/2E. Although the potentials Va are typically very small, so are Am'^/2E; if Va are 
comparable to or larger than Am'^/2E, matter can strongly affect neutrino oscillations. 



8.1 Evolution equation 

We shall now consider neutrino oscillations in matter in some detail. Neutrinos of all three 
flavours - z/g, z/^ and u^- - interact with the electrons, protons and neutrons of matter 
through neutral current (NC) interaction mediated by bosons. Electron neutrinos in 
addition have charged current (CC) interactions with the electrons of the medium, which 
are mediated by the exchange (see fig. 1). 
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Figure 5: Neutrino scattering diagrams 



Let us consider the CC interactions. At low neutrino energies, they are described by 
the effective Hamiltonian 

HCC = ^ - 75)i^e][^e7'^(l " 75)e] = ^ [67^(1 - 75)e][/>e7^(l - 75)^e] , (84) 

where we have used the Fierz transformation. In order to obtain the coherent forward 
scattering contribution to the energy of z/g in matter (i.e. the matter-induced potential 
for Ug) we fix the variables corresponding to z/g and integrate over all the variables that 
correspond to the electron: 

Hejfiye) = {Hcc) electron = ^'eK^^e • (85) 

Furthermore, we have 

(e7oe) = (e^e) = Ne , (e7e) = (v^) , (e7o75e) = {^^) , (e775e) = (o-g) , (86) 

J^e 

where A^e is the electron number density. For unpolarized medium of zero total momentum 
only the first term survives, and we obtain 

{y,)cc = Vcc = V2GFN,. (87) 

Analogously, one can find the NC contributions Vnc to the matter-induced neutrino poten- 
tials. Since NC interaction are flavour independent, these contributions are the same for 
neutrinos of all three flavours. In an electrically neutral medium, the number densities of 
protons and electrons coincide, and the corresponding contributions to Vnc cancel. The di- 
rect calculation of the contribution due to the NC scattering of neutrinos off neutrons gives 
{Va)NC = —GpNn/V^, where is the neutron number density. Together with eq. ( |57| ) 
this gives 

K = v^G^ (A^e-^) , V, = Vr = V2Gp [-^y (88) 

For antineutrinos, one has to replace Va —Va- 

Let us now consider the evolution of a system of oscillating neutrinos in matter. In 
vacuum, the evolution is most easily followed in the mass eigenstate basis. In matter it is 
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more convenient to do that in the flavour basis because the effective potentials of neutrinos 
are diagonal in this basis. Consider first the two-fiavour case. As usual, we write Ufi = Uv^ 
where Vfi and Vm are two- component vectors of neutrino fields in the fiavour and mass eigen- 
state bases and the matrix U is given by eq. (0). In the absence of matter, the evolution 
equation in the mass eigenstate basis is i{d/dt)\h'rn) = Hm\i^m), where Hm = diag{Ei, E2). 
This gives the evolution equation in the fiavour basis: i{d / dt)\v fi) = Hfi\ufi) = UHmU^Ufi) . 
For relativistic neutrinos Ei c:^ p + mf/2E, and we thus obtain 




{p+^)-^cos2eo ^ sin 2^0 

^ sin 2^0 (p + %^) + ^cos2^o 




^9) 




Here and stand for time dependent amplitudes of finding the electron and muon neu- 
trino respectively. The expressions in the brackets in the diagonal elements of the effective 
Hamiltonian in eq. ( p9D coincide. They can only modify the common phase of the neu- 
trino states and therefore have no effect on neutrino oscillations which depend on the phase 
differences. For this reason one can omit these terms. The evolution equation describing 
neutrino oscillations in vacuum in the fiavour basis then takes the form 

cos 2^0 ^sin2^^o\/^e\ .gp. 
^ sin 2^0 ^cos2eJ[uJ- ^''^ 

We now proceed to derive the neutrino evolution equation in matter. To do that, one 
has to add the matter-induced potentials VJ= and to the diagonal elements of the effective 
Hamiltonian Hfi in eq. (pOl). Notice that Vg and contain a common term due to NC 
interactions. As we already know, such common terms in the diagonal elements are of no 
consequence for neutrino oscillations; we can therefore omit them Q This gives 

■±( ^e\^( -^cos2e, + V2GfN, ^sin2^^o \( .qn 
dt\v,) \ ^ sin2eo ^ cos2^^o / U J ' 

This is the evolution equation which describes z/g ^ z/^ oscillations in matter. The equation 
for Vf. ^ oscillations has the same form. In the two-flavour approximation, ^ Vj. 
oscillations are not modifled in matter since = Vr] however, in the full 3-flavour framework 
matter does influence the t'^ Ur oscillations because of the mixing with u^, see sec. |0 . 



8.2 Constant density case 

Let us now consider the evolution equation (0). In general, the electron number density Ne 
depends on the coordinate along the neutrino trajectory or, in our description in eq. (0), on 

^The NC contribution does, however, affect the oscillations between the usual and "sterile" (electroweak 
singlet) neutrinos in matter. 
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t. We shall first consider a simple case of constant matter density and chemical composition, 
i.e. A'^e = const. Diagonalization of the effective Hamiltonian in ( PT| ) gives the following 
neutrino eigenstates in matter: 

Pa = J^e cos 9 + Ufj, sin ^ , 

ub = —Uf. sinO + cos9 , (92) 



where the mixing angle 6 is given by 



sin 26 Q 



tan 29 = -j—^ ^ . (93) 

^cos2eo- V^G^iVe 

It is different from the vacuum mixing angle and therefore the matter eigenstates and 
vb do not coincide with mass eigenstates vi and z/2 . The difference of neutrino eigenenergies 
in matter is 



Ea-E 



B 



where 



\ \~2F --^Gf NA + i-l^] sin' 2^0 • (94) 

It is now easy to find the probability of z/g ^ oscillations in matter: 

P{ve v^- L) = sin' 29 sin' [vr-^] , (95) 

27r 2n 

~ y cos 2^0 -V2GfN:)\ (^)%in' 29^ 

It has exactly the same form as the probability of oscillations in vacuum (^), except that 
the vacuum mixing angle 9q and oscillation length Ig^c are replaced by those in matter, 9 
and Im- In the limit of zero matter density 9 = 9^, Im = loso and the vacuum oscillation 
probability is recovered. 

The oscillation amplitude 

(^)%in'2eo 

sin' 29 = 2 2 (97) 

(^cos2^o- v^Gi.iV,) +(^) sin' 2^0 

has a typical resonance form, with the maximum value sin' 29 = 1 achieved when the 
condition 

V2GfN, = ——cos29o (98) 

is satisfied. It is called the MSW resonance condition. From (|93| ) or (|97|) it follows that 
when (^8]) is fulfilled, mixing in matter is maximal {9 = 45°), independently from the vacuum 
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mixing angle ^o- Thus, the probabihty of neutrino flavour transition in matter can be large 
even if the vacuum mixing angle is very small! 

For the resonance enhancement of neutrino oscillations in matter to be possible, the 
r.h.s. of (p8| ) must be positive: 

Am^ cos 260 = {ml - m?) (cos^ 60 - sin^ ^0) > , (99) 

i.e. if 1^2 is heavier than z/i, one needs cos^^o > sin^6'o, and vice versa. It follows from 
eq. that the condition (^) is equivalent to the requirement that of the two mass 

eigenstates ui and z/2, the lower-mass one have a larger z/g component. If one chooses the 
convention cos26'o > 0, (as is usually done) then (|99|) reduces to Am^ = Aml^ > 0. The 
resonance condition for antineutrinos is then Am^ < 0. Therefore, for a given sign of Am^, 
either neutrinos or antineutrinos (but not both) can experience the resonantly enhanced 
oscillations in matter. 



8.3 Adiabatic approximation 

Let us now discuss the realistic case of matter of varying density. Typically, one deals with 
situations when a beam of non-monochromatic neutrinos (i.e. neutrinos with some energy 
distribution) propagates in a medium with certain density proflle. If Am^ is of the right 
order of magnitude, then for any value of the matter density (or a signiflcant portion of 



them) there is a value of neutrino energy for which the resonance condition (98) is satisfled. 
Conversely, every value (or a significant portion of the values) of neutrino energy E "finds" a 
value of the matter density for which the resonance condition ( pHD is satisfied. If the neutrino 
beam is monochromatic (as, e.g., in the case of the solar ^Be neutrinos), the resonance 
enhancement of neutrino oscillations is still possible if the corresponding resonance density 
is within the density range of the matter in which neutrinos propagate. Thus, the MSW 
resonance condition does not involve any fine tuning. 

In general, for oscillations in a matter of an arbitrary non-uniform density, the evolution 
equation (|9l|) does not allow an analytic solution and has to be solved numerically. However, 
there is an important particular case in which one can get an illuminating approximate 
analytic solution. This is the case of slowly (adiabatically) varying matter density. 

We shall start with a semi-quantitative discussion. Consider electron neutrinos born 
in a matter of a very high density, far above the MSW resonance one (e.g., in the core 
of the sun). We shall assume that the matter density is monotonically decreasing along 
the neutrino trajectory. From eq. ( pSf ) it follows that the mixing angle in matter at the 
neutrino production point 6i ~ 90°, which means that the neutrino mixing is strongly 
suppressed by matter. As neutrinos propagate towards the regions of smaller density, their 
mixing increases (mixing angle 6 decreases); it becomes maximal at the resonance point, 
where 6 = 45°. As neutrinos propagate further towards smaller densities, their mixing 
angle continues decreasing; it reaches the value 9f = 9q at densities <C {Ne)Msw, where 
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Figure 6: Neutrino energy levels in matter vs electron number density Ng. Dashed line - in the 
absence of mixing, solid line - with mixing. 

{Ng)Msw is the resonance value of the electron number density given by eq. (pS]). From 
eq. ( P^ it follows that at the production point, where 9 ~ 90°, the produced almost 
coincide with the matter eigenstate neutrinos z/^. If the matter density changes slowly 
enough (adiabatically) along the neutrino path, the neutrino system has enough time to 
"adjust" itself to the changing external conditions. In this case the transitions between the 
matter eigenstates ua and ub are exponentially suppressed. This means that if the system 
was initially in the eigenstate ub, it will remain in the same eigenstate. However, the flavour 
composition of this eigenstate changes as neutrino propagates in matter because the mixing 
angle 6 that determines this composition is the function of matter density. Since at the final 
point of neutrino evolution 9 9q, eq. (^) tells us that the matter eigenstate neutrino ub 
at this point has the component of originally produced z/g with the weight sin^ 9q and the 
component of z/^ with the weight cos^ i-e. the transition probability is 

P{i^e = cos^ ^0 . (100) 

This means that in the case of small vacuum mixing angle, one can have almost complete 
adiabatic conversion of z/g to z/^! 

This is illustrated by fig. ^ which shows the energy levels of and ub along with those 
in the absence of mixing (i.e. of z/g and u^) as the function of the electron number density. 
In the absence of mixing the energy levels cross at the MSW resonance point, but with 
nonvanishing mixing the levels "repel" each other, and the avoided level crossing results. 
If the probability of the transition between the two matter eigenstates is small, neutrinos 
produced as z/g at high densities and propagating towards smaller densities follow the upper 
(ub) branch and end up on the level that corresponds to z/^ at small Ne- This resonance 
conversion is similar to the well known Landau - Zener phenomenon in atomic and molecular 
physics. 

The expression ( p,00| ) for the conversion probability looks paradoxical: the smaller the 
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vacuum mixing angle, the larger the probability that the initially produced Ve will be con- 
verted into or vice versa. Does this mean that in the limit of vanishing Qq one can still 
have strong neutrino conversion? The answer, of course, is no. The reason is simple: if 
6*0 becomes too small, the adiabaticity of the conversion gets broken, and eq. ( |100| ) ceases 
to be valid. Similar situation takes place in the case of neutrino oscillations in a matter 
of constant density: if the MSW resonance condition is satisfied, the oscillation amplitude 
sin^ 26'o = 1, no matter how small ^o! however, in the limit ^ the phase of the second 
sin^ factor in eq. (|95|) vanishes, and no oscillations occur. 



We shall now turn to a more quantitative description of the neutrino conversion in the 
adiabatic regime, which will also allow us to establish the domain of applicability of the 
adiabatic approximation. At any instant of time, the effective Hamiltonian Hfi{t), which is 
given by the matrix on the r.h.s. of eq. (|9T|), can be diagonalized by a unitary transformation 

Ufi = U{t)u , U{tyHfi{t)U{t) = Hd{t) = dtag{EA{t), Esit)) . (101) 

Here u = (z^a(^), ^B(t))'^ is the vector of the instantaneous matter eigenstates, -Ea(^) and 
Esit) are instantaneous eigenvalues of Hfi{t), and the matrix U{t) has the form (|56|) except 
that the vacuum mixing angle 6*0 has to be replaced by the mixing angle 9(t) given by eq. ( p3| ) 
with Ne = Nf,(t). The evolution equation in the basis of the instantaneous eigenstates can 
therefore be written as i{d/dt)v = [Hd — iW{dU /dt)]^, or 




EA{t) -id{t) 




(102) 



where 6 = dO/dt. Notice that the effective Hamiltonian in this basis is not diagonal since the 
mixing angle 9 is not constant, i.e. the matter eigenstate basis changes with time. If the off- 
diagonal terms are small, i.e. \9\ <^ \Ea — Eb\, the transitions between the instantaneous 
eigenstates z/^ and ub are suppressed. This corresponds to the adiabatic approximation. 
The adiabaticity condition can be written as 



2\e\ sin 2^0 



-'-'-n^-^:zi&l''c.cl«i. (103) 



where Ea — E^ and Vcc are given by eqs. (p^ and (^7]) respectively. The parameter 7 is 
called the adiabaticity parameter. In the adiabatic limit, the effective Hamiltonian in ( |102| ) 
is diagonal, and time evolution of the matter eigenstates is very simple - they just receive 
phase factors. 

Suppose that at a time t = ti the electron neutrino was born: 

z/(tj) = z/g = cos 6*4 va + sin^^j ub ■ (104) 
Then, in the adiabatic approximation, at a time we have 

/, \ /I -i EA(t')dt' . „ -i f^f EB(t')dt' 

u{tf) = cosOie ■''i ^' UA + smOie ■''^ ^' vb ■ (105) 
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Taking into account that aXt = tf the mixing angle 9{tf) = Of is different from 9i, one then 
finds 

P{ve ^ z^/J = ^ - 2 cos cos 26 f - - sin 2^^ sin 20/ cos $ , (106) 

where 

$= f\EA-EB)dt'. (107) 

The second term in eq. ( |106| ) is a smooth function of tj, while the third term oscillates 
with time. If the matter density at the neutrino production point is far above the MSW 
resonance one, then sin20j ^ and the third term is strongly suppressed. In this case the 
non-oscillatory neutrino conversion takes place, with the probability P{ve v^) = cos^0/. 
If, in addition, the final point of neutrino evolution corresponds to a very low matter density. 
Of = 9o and eq. ( |100| ) is recovered. Notice that eq. ( |106| ) is symmetric with respect to the 
interchange of the initial and final evolution points ti ^ tf. 

One can readily generalize this analysis to include possible transitions between the mat- 
ter eigenstates pa and pb due to the violation of the adiabaticity. Let P' be the probability 
that the va ^ transitions have occured in the course of the evolution of the neutrino 
system ( "hopping probability" ) . Then 



P{v, ^v^)^--- co^lQi cos2%(l - 2P') , (108) 

where we have omitted oscillating terms which average to zero. In the case of very small 
vacuum mixing angle and (A^e)/ < {'Ne)msw < (A^e)i (or (A^e)j < {J^e)msw < (^e)/), it 
reduces to a very simple expression 

P(z/e^Z/^)~l-P'. (109) 

In the Landau-Zener approximation the hopping probability P' can be written as 

where 7^ is the adiabaticity parameter taken at the MSW resonance point. In the adiabatic 
limit, 7r ^ 1, one has P' ~ and eq. ( |106| ) is recovered. In the non-adiabatic limit 7^ -C 1, 
P' ~ 1, and the survival and transition probabilities interchange; in that case the transition 
probability is small when is small. It should be noted that in the extreme non-adiabatic 
limit eq. ( |110| ) for the hopping probability is not quite correct and more accurate expressions 
should be used. It is important that the transition probability can be quite sizable even if 
the adiabaticity parameter is not very large: for 7^ ~ 1 it can be as large as about 80%, 
and one can, e.g., consider the condition 7^ > 3 (which can result in as large transition 
probabilities as 99%) as the adiabaticity criterion. 

Let us now discuss the adiabaticity condition ( |103| ) in more detail. Since the difference 
of the eigenenergies reaches its minimum equal to {Am^/2E) sin 26*0 at the MSW resonance 
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(see eqs. (0) and (p8|) and fig. H), the adiabaticity is worst at the resonance. If the 
adiabaticity condition at the resonance is satisfied, it will then also be satisfied away from 
it []. The adiabaticity parameter at the MSW resonance is 



/ Am^ ^ \ ^ 1 sin^ 26*0 Am^ 
7r = sin 2^0 -■ = Lp, 111 

^ V V \Vcc\res cos 2^0 2E ' ' ^ ' 



where Lp = \Vcc/ycc\res = l-^e/^^el^es the density scale height at the resonance, i.e. the 
characteristic distance over which the electron number density varies significantly in the 
resonance region. 

The adiabaticity condition has a simple physical meaning. Let us define the resonance 
width at half height Ar as the spatial width of the region where the amplitude of neutrino 
oscillations in matter sin^ 2^ > 1/2. It is easy to obtain 

Ar~2tan2^o^P (112) 

From ( P^l ) we find that the oscillation length at the resonance is {lm)res = '^'^/lEA — Eslres = 
{AtiE/ Am?)/ sin 2^0- Therefore the adiabaticity parameter at the resonance ( |1 1 1| ) can be 
written as 

Ar 

7r = TT , (113) 

\''m)res 

i.e. the adiabaticity condition 7^ > 3 is just the condition that at least one oscillation length 
fit into the resonance region. 

As follows from (^) and ( |1 1 1| ) , both the MSW resonance condition and the adiabaticity 
parameter at the resonance 7^ depend on neutrino energy. Therefore the efficiency of the 
matter-enhanced neutrino fiavour conversion is energy dependent. This should have impor- 
tant consequences, for example, for solar or supernova neutrinos, leading to characteristic 
distortions of the energy spectra of detected neutrinos. 

For neutrino oscillations in a matter of density which decreases monotonically along 
the neutrino path (e.g., for neutrino oscillations in the sun), the average neutrino survival 
probability Ps is schematically shown in fig. (^. It has a typical "bathtub" shape. In 
the case of small vacuum mixing angle, the low energy edge of the curve corresponds to 
no suppression (resonance density exceeds the maximal available matter density {Ne)max ^ 
there is no MSW resonance). The steep decrease of Ps corresponds to the value of E / Am? 
at which the MSW resonance condition ( pS] ) is satisfied at N^. = {Ne)max- The bottom of the 
suppression bathtub corresponds to the adiabatic regime, whereas the slow increase of Ps 
at higher energies is due to the loss of adiabaticity. The high energy edge of the suppression 
curve corresponds to the non-adiabatic regime. In the case of larger vacuum mixing angles, 
the suppression curve shows a similar behavior, but with smaller amplitudes of changes of 
the survival probability Ps and a wider adiabatic region. 

^Assuming that |^| does not have strong maxima outside the resonance region. 
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Figure 7: Averaged survival probability vs Ej in the case of neutrino oscillations in a matter 
of monotonically decreasing density. Solid line - small dashed line - large ^o- 



It should be noted that the fulfillment of the adiabaticity condition ( [103| ) is not, in 
general, sufficient for the strong flavour conversion to take place: in addition, one has to 
require that the total width (column density) d of the matter traversed by neutrinos be 
large enough. If matter density changes slowly along the neutrino trajectory, the scale 
height Lp is large, and so the adiabaticity parameter 7^ and the resonance width Ar given 



by eqs. ( |111|) and (|112|) can formally be very large, too; however, in reality the resonance 
width cannot exceed the total width of the matter slab traversed by neutrinos. For the 
transition probability to be appreciable, the column density d must satisfy 



d 



dt' > 



Gi? tan 26*0 



(114) 



Problem 16. Derive eq. l \106i ). 

Problem 17. Derive eq. dllSj ) for the MSW resonance width. 

Problem 18. Assuming that the ^ Vb transitions occur (with probability P') in a narrow 
region around the MSW resonance and that the change of the oscillation phases in that 
region is small, derive the non-adiabatic analog of eq. l \10(\ ). Show that upon averaging of 



the oscillating terms eq. ( \10^ ) is obtained. 
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8.4 3 flavour case 



The evolution equation describing neutrino oscillations in matter in the case of three neutrino 
flavours i/g, i'^ and Vr is 



d 
dt 
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(115) 



Here U is the 3-fiavour vacuum mixing matrix defined in eq. (0) or (^), and Vcc was 
defined in eq. (^). In general, it is rather difficult to study this equation analytically, even 
in the adiabatic approximation. There are, however, several important particular cases when 
the evolution of the 3-flavour neutrino system can be approximately expressed through the 
2-flavour amplitudes. One such case is the case of the hierarchical mass squared differences 
([66|). It is of particular interest because this is exactly the type of the hierarchy which is 
suggested by the solar and atmospheric neutrino experiments. We shall now discuss this 
case in some detail. 

The mixing matrix U in the three-flavour case can in general be written as (see, e.g., 
) 



U = V2^WnVi2D = VD . 



(116) 



where 
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and D = diag{e 1, e *'''^). The phase S in Wis is the usual Dirac-type CP violating 
phase, whereas the phases ipi and ip2 are present only in the Majorana case. It immediately 
follows from eq. (|115|) that the Majorana phases have no effect on neutrino oscillations, and 
therefore one can omit the factor D and write U = V. In this case the matrix U defined in 
( |116| ) coincides with the one in (p5|). 

We shall first consider the case (|67D which is relevant for atmospheric, reactor and 
accelerator neutrino experiments. Subtracting the term from all the diagonal elements 
of the matrix = diag{m\, m\, m|) and neglecting Am^i compared to Am|]^ we get — > 

= diag{0, 0, Am|^). The new matrix commutes with the rotation matrix Vu, so 
the evolution equation ( |115| ) can be rewritten as i{d/dt)i' = [V23Wi3{M'^ /2E)WisV23 + V]^, 
where V is the second matrix on the r.h.s. of eq. ( p.l5| ): V = diagiVcc, 0, 0). Since V 
commutes with V23, it is convenient to go into a new basis defined through u = V23I'. The 
evolution equation in this basis is 
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i.e. the problem is reduced to a 2-flavour problem of oscillations between ui = and 
= -523 i^^i + C23 I'r- These oscillations are governed by the large mass squared difference 
Amli and mixing angle ^13. The CP- violating phase 6 is unobservable in the limit Am^i = 
and so we have omitted it. It is convenient to subtract the term {Amli/2E + Vcc)/'^ 



from all the diagonal elements of the effective Hamiltonian in eq. ( |118| ), so as to make the 
Hamiltonian of the two-flavour system (z/i, z/3) traceless. The two-flavour evolution equation 
then takes the form 




cos 2^,3 + ^ ^ sin 2^,3 
^sin2^,3 ^cos2^^,3-^ )[-^)' ^ ^ 



For an arbitrary initial state ?'(to) the solution at a time t can be written as z/(t) = 
5'(t, to)^'(O), where S(t,to) is the evolution matrix that corresponds to the effective Hamil- 
tonian in eq. (|119|) . The unitary 2x2 matrix S'(t, to) can be parametrized as 



S{t,to) = Z -irW , (120) 

where r are the Pauli matrices in the (i>i, z/3) space, and the real parameters Z and W = 
(Wi, W2, W3) satisfy Z'^ + 'W'^ = 1. The matrix S(t, to) can be found for any matter density 
proflle using the standard techniques developed for the 2-flavour neutrino oscillations. In 
particular, for long-baseline accelerator experiments, neutrino path goes through the mantle 
of the earth where to a very good approximation the matter density can be considered to be 
constant. In this case the parameters Z and W can be trivially found. For the atmospheric 
neutrino oscillations, in the case of neutrinos traversing both the mantle and the core of 
the earth, simple analytic expressions for the parameters Z and W can be obtained in the 



step-function approximation to the density proflle of the earth |^ . With known Z and W 
it is easy to calculate the probabilities of neutrino oscillations: 

Pil^e ^ l^f.) = S^3 P2 , P{l^e ^ Vr) = 4 ^2 , (121) 

P{y^,^yr) = s^3C^3[2-P2-2(Zcos<l> + lV3sin$)], (122) 

where 

^ = \f{^ + ycc\dt', P, = W^ + W^, (123) 



2 Jto \ 2E 



P2 being the probability of the ui transitions described by eq. (|119|) . It is not difficult 
to take into account corrections due to nonvanishing Am'^i ||3^ - one just has to add the 



term {—Aml^/E) to the integrand in eq. (|123|) . Eqs. (|121|) - ( p.23| ) generalize the results 



([69| ) - ( [71]) to the case of neutrino oscillations in matter. 

Problem 19. Derive eqs. ( \121\ ) - (\12^ ). For the case of neutrino oscillations in a matter of 
constant density, find explicit formulas for the parameters Z and W. Show that in the limit 
of vanishing matter density eqs. - ( |yi|^ are recovered. 
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In the limiting case (^) which is relevant for the solar neutrino oscillations, and assum- 
ing, as before, the hierarchy (B^), one finds that the z/g survival probability is again given by 



eq. (|7J). However in the presence of matter the two- flavour survival probability P in this 
expression is no longer given by eq. ([75|); it has to be calculated from the 2-flavour evolution 
equation in matter of the type (|9T|) governed by the mass squared difference Am^ = Am^i 
and mixing angle 9q = 9i2, except that the matter- induced potential for neutrinos Vcc must 
be replaced by the effective one Veff = cf^ Vcc @- 

Problem 20. Assuming that lS.m\j2E ~ I^ml^l^E > I\m\j2E ~ Vcc, {l^m\j2E)L > 1, 
derive eq. ( fT^ j in the case of neutrino oscillations in matter and show that the survival 
probability P in it has to be found from the 2-flavour evolution equation with Veff = 0^3 Vcc- 
Hint: consider neutrino evolution in the basis u' = Vi^V2^i' and take into account that the 
evolution of the state z/3 decouples. 

Problem 21. Derive expressions for 3-fiavour oscillation probabilities in terms of the 
relevant 2-flavour evolution matrix in the case Ue3 = 0. Find explicit formulas for the case 
of the oscillations in a matter of constant density. Show that in the limit of vanishing density 
eqs. (f7^)- ^7^ ) are recovered. 



As we have already discussed, matter affects neutrino and antineutrino oscillations dif- 
ferently: for cos26'o > and Am^ > 0, it can enhance oscillations of neutrinos but always 
suppresses oscillations of antineutrinos, whereas in the case Am^ < the situation is op- 
posite. This means that even in the absence of the CP-violating phase 6 in the mixing 
matrix (|65|) , the probabilities of neutrino oscillations are different from those of antineu- 
trino oscillations, i.e. matter induces the CP violating effects. This is because the usual 
matter is itself CP (and also C and CPT) asymmetric since it consists of particles and not 
of antiparticles, or in general not of equal number of particles and antiparticles. This CP 
violation is, however, "macroscopic" rather than microscopic; although it may be useful for 
experimental searches of matter effects on neutrino oscillations, it does not say anything 
about the CP violation in the leptonic sector at the fundamental level. Thus, matter ef- 
fects can mimic the true CP violation and therefore may make it difficult to disentangle 
the genuine CP violation from the macroscopic one in the long baseline neutrino oscillation 



experiments [39 



There is, however, a circumstance that may in principle facilitate experimental searches 
of the genuine CP violation in the leptonic sector. In vacuum, CP-odd and T-odd asym- 
metries of the neutrino oscillation probabilities coincide; they both depend on the same 
CP- and T-violating parameter 6. For neutrino oscillations in a medium, CP is violated by 
the matter but T in general is not: it is only violated if the matter density profile is T- 
asymmetric, i.e. asymmetric with respect to the interchange of the positions of the neutrino 
source and detector. A matter with T-symmetric density profile (and in particular, a matter 
of constant density) does not induce any T violation, and an experimental observation of a 
T violating asymmetry P(i^a ^b) -~ P{^b ^a) in such a matter would directly measure 
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the fundamental T (and CP) violating phase 5. 



Matter effects on propagating neutrinos are not limited to the MSW effect discussed in 
this section. Matter can also modify significantly neutrino decays in the case of unstable 
neutrinos and neutrino spin and spin-flavour precession in magnetic fields In par- 



ticular, the spin-flavour precession (simultaneous change of neutrino chirality and flavour 
due to interaction of neutrino flavour-off-diagonal magnetic moments with magnetic fields) 
can be resonantly enhanced in matter; this phenomenon is similar to the MSW effect. Neu- 
trino oscillations in matter can also be amplified even if the neutrino mixing angles both 
in vacuum and in matter are small, i.e. if there is no MSW effect (no level crossing). This 
is possible if the change of matter density along the neutrino path and the change of the 

I7|. In this case the oscilla- 

This 



oscillation phase are "synchronized" in a certain way |]42|, H3, 



tion probability can be strongly enhanced due to a cumulative build-up effect ||42| , ^ 



phenomenon is of the same nature as the parametric resonance in mechanical or electro- 
magnetic oscillations. It has been discovered recently that this parametric enhancement of 
neutrino oscillations can lead to observable effects in the oscillations of atmospheric ^ 
and solar iHSl E6l neutrinos in the earth. 



9 Double beta decay 

In some cases when the ordinary beta decay processes (|l|) are energetically forbidden, the 
double beta decay processes, in which a nucleus A{Z, N) is converted into an isobar with 
the electric charge differing by two units, may be allowed: 

A{Z, N) A{Z ± 2, ^ 2) + 2e^ + 2z/e(2z/e) . (124) 

In such decays two neutrons of the nucleus are simultaneously converted into two protons, 
or vice versa. At the fundamental (quark) level, these are transitions of two d quarks into 
two u quarks or vice versa (see fig. ^ja). Double beta decays are processes of the second 
order in weak interactions, and the corresponding decay rates are very low: typical lifetimes 
of the nuclei with respect to the 2(3 decay are T > 10^^ years. The processes ( p.24| ) are 
called 2/32i/ decays. Two-neutrino double beta decays with the emission of two electrons 
(2/3~) were experimentally observed for a number of isotopes with the half-lives in the range 
10^^ — 10^^ years; there are very few candidate nuclei for 2/3"'" decay, and the experimental 
observation of this process is difficult because of the very small energy release {Q values). 

If the lepton number is not conserved, the electron neutrino or antineutrino emitted in 
one of the elementary beta decay processes forming the 2/3 decay can be absorbed in another 
(fig. Pb), leading to the neutrinoless double beta decay (2/?0z/): 

A{Z,N) ^ A{Z ±2,NT'2) + 2e^ . (125) 

Such processes would have a very clear experimental signature: since the recoil energy of a 
daughter nucleus is negligibly small, the sum of the energies of the two electrons or positrons 
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in the final state should be equal to the total energy release, i.e. should be represented by 
a discrete energy line. Therefore 2/30z/ decays could serve as a sensitive probe of the lepton 
number violation. In some extended models exotic modes of 2/30i^ decay are possible, e.g. 
decays with Majoron emission. In this case the sum of the energies of two electrons or 
positrons is not a discrete line, but their spectrum is expected to be different from that with 



no Majoron emission ^6 
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Figure 8: 2/3 decay 



Neutrinoless 2(3 decays break not only the lepton number; since the absorbed z/g or 
Vf. has a "wrong" chirality, 2/?0z/ decays also break chirality conservation (this is another 
example of the chiral prohibition rule). Therefore 2/50z/ decay is only possible if neutrinos 
have nonzero Majorana mass. Indeed, if 2/30z/ decay is mediated by the standard weak 
interactions, the amphtude of the process is proportional to neutrino mass (fig. |b): 

A(2/50z/) oc Ul = (m.Je// • (126) 

i 

Notice that this expression contains U^^ rather than |t/eiP- If CP is conserved in the leptonic 
sector, the mixing matrix U^^ can always be made real; however in this case the mass 
parameters in (|126|) (eigenvalues of the neutrino mass matrix) can be of either sign, 
their relative signs being related to the relative CP parities of neutrinos. This means that 
in general significant cancellations between various contributions to the sum in ( |1 26| ) are 
possible. As we discussed in sec. ^, a pair of Majorana neutrinos with equal physical masses 
|mi|, opposite CP parities and maximal mixing is equivalent to a Dirac neutrino. It is easy 
to see that such a pair does not contribute to the amplitude in ( |126|) . Analogously, the 
contribution of a quasi-Dirac neutrino is strongly suppressed. 

Naively, one could expect that the rates of 2/?0z/ decays to be much higher than those 
of 2j32v decays because the final state in 2/30z/ decays contains fewer particles than that 
of 2j32v decays. In reality the situation is more complicated - the rates of 2/?0z/ decay 
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are strongly suppressed by a very small factor {my^)1jj/El, where £"0 (~ a few MeV) is 
the typical energy release in 2[3 decay. Neutrinoless double beta decay was searched for 
experimentally but up to now have not been discovered. The experiments allowed to put 
the upper bound on the effective Majorana neutrino mass, the best limit coming from the 
Heidelberg - Moscow experiment on 2(3 decay of ^^Ge p6l: 



(m,Je// < 0.2-0.6 eV, (127) 

depending on the value of the nuclear matrix element which is not precisely known. If the 
2j3Qv decay is discovered, it will be possible to infer the value of the effective Majorana 
neutrino mass (w^i/e)e//- As follows from ( |126| ), this would give the lower limit on the mass 
of the heaviest neutrino. 

In extended gauge models such as the left right symmetric (LR) or SUSY models, 
additional mechanisms of 2/30z/ decay are possible. For the LR models a typical diagram 
is shown in fig. ||c. It appears that no Majorana mass rriL of ul is necessary in this case, 
i.e. 2/?0z/ decay can occur even if rriL = 0. This, however, is incorrect: in all extended 
models in which 2f30i' decay occurs, the Majorana masses of are different from zero. In 
particular, in the LR models in which the diagram ||c exists, the seesaw mechanism operates 
and so ttil 7^ 0. It can be shown in a model-independent way that 2j30u decay always means 
rriL 7^0 



10 The solar neutrino problem 

It is generally believed that the stellar energy is provided by thermonuclear reactions. For 
so called main sequence stars (to which our sun belongs), the main result of these reactions 
is a fusion of hydrogen into helium: 

4p + 2e" = ^He + 2i/e + 26.73 MeV . (128) 

The main reaction chain in which the fusion process ( |128| ) occurs in the sun is the pp cycle 
(fig. whereas the carbon-nitrogen-oxygen (CNO) cycle is responsible for less than 2% of 
solar energy. Neutrinos are emitted in six reactions of the solar pp cycle (Table |ip. They 
are produced in nuclear beta decay or electron capture reactions, and their energy spectra 
are well known (see fig. |l^). Three of these reactions - pep and the electron capture into 
two different final states of ^Be - produce monochromatic lines, whereas the neutrinos 
born in the other three reactions have continuous energy spectra. The fluxes of the solar 
z/g are calculated in the framework of the standard solar models. These models are based 
on the assumptions of local hydrostatic equilibrium, thermonuclear nature of solar energy, 
energy transfer in the sun by radiation and convection and the requirement that the present 
values of the solar radius, mass, luminosity and He/H ratio be reproduced as a result of the 
solar evolution. There are about 20 different solar models developed during the last decade 
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by more than ten different groups of authors, all based on the above mentioned principles. 
All these models except one agree with each other well (see fig. [TT]). A number of recent 
models also took into account diffusion of helium and heavy elements. These models were 
able to reproduce with a very high accuracy the solar sound velocities inferred from the 



helioseismological measurements [Q. This gives further credibility to the standard solar 
models. 

The calculated fluxes of different components of the solar v^, flux depend not only on 
the general principles laid in the basis of the standard solar models, but also on a number 
of input parameters, such as the solar opacity and various nuclear cross sections. These 
parameters are known to a varying degree of accuracy . The flux of the most abundant 
pp neutrinos is closely related to the well known luminosity of the sun, and therefore has a 
very small uncertainty, ~ 1%. On the other hand, the flux of the energetic ®B neutrinos has a 
relatively large uncertainty, ~ 20%. The flux of the most energetic hep neutrinos is the least 
certain one - it is only known to an accuracy of 1 - 2 orders of magnitude. Fortunately, this 
flux is very small, and its contribution to the detection rates of solar neutrinos is practically 
negligible. It can, however, affect the high energy part of the solar neutrino spectrum. 

The solar neutrinos have been detected in flve solar neutrino experiments. Historically, 
the flrst one was the Homestake experiment of Davis and his collaborators. It is based on 
the reaction 

i/e + ^^Cl ^ ^^Ar + e" . (129) 

The radioactive argon atoms produced in this reaction are extracted using chemical methods 
and counted in proportional counters. The energy threshold of reaction ( |129D is 0.814 MeV, 



so only the and ^Be and pep neutrinos are detected in the Homestake experiment, the 



largest contribution coming from the B neutrinos (see flgs. |T0| and |T^). Radiochemical 
techniques are also used in two other solar neutrino experiments - SAGE and Gallex, in 
which the reaction 

z/e + ^^Ga^^^Ge + e". (130) 
is employed. The energy threshold of this reaction is 0.234 MeV, and so the gallium exper- 



Table 1: Solar neutrino production in the pp chains (from ref. |[l3[|). 
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Figure 9: Energy generation in the Sun via the pp chains (from ref. 




Figure 10: Solar neutrino spectrum and thresholds of solar neutrino experiments as indicated 
above the figure (from ref. ||4^ ). 
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iments can also detect the lowest energy pp neutrinos. Since the flux of the pp neutrinos 
is very large, they are expected to give the main contribution to the event rates in the 
SAGE and Gallex detectors (flg. |12D. The remaining two experiments - Kamiokande and 
its up-scaled version Super-Kamiokande - are the water Cherenkov detectors and use the 
neutrino-electron scattering reaction 

i^a + Ua + e~ . (131) 

to detect solar neutrinos. This reaction has zero physical threshold, but one has to introduce 
energy cuts to suppress the background. In the Kamiokande experiment solar neutrinos with 
the energies E > 7.5 MeV were detected, whereas the threshold used by Super-Kamiokande 
is at present 5.5 MeV. With these energy cuts, the Kamiokande and Super-Kamiokande 
detection rates are only sensitive to the component of the solar neutrino flux. Reaction 
(|131|) has a very interesting feature - for neutrino energies E ^ rrie the angular distribution 
of the recoil electrons is forward-peaked, i.e. it points in the direction of the momentum of 
the incoming neutrino. The angular distributions of neutrinos detected in the Kamiokande 
and Super-Kamiokande experiments have a prominent peak at 180° from the direction to 
the sun, which is a beautiful proof of the solar origin of these neutrinos. The origin of the 
neutrinos detected in the radiochemical experiments is unknown, and our belief that they 
come from the sun is based on the fact that no other sources of sufficient strength and 
energy are known. 

In all five solar neutrino experiments, fewer neutrinos than expected were detected, the 
degree of deflciency being different in the experiments of different types (flg. |1^). This is 
the essence of the solar neutrino problem. The observed deflciency of solar neutrinos could 
be due to an insufficient knowledge of solar physics or an error in some input parameters 
(astrophysical solution), or due to experiment-related errors (such as miscalculated detection 
efficiency or cross section), or due to some unknown neutrino physics (particle physics 
solution of the solar neutrino problem). As we have already discussed, the standard solar 
models are well grounded, agree with each other and with helioseismological observations. 
All the solar neutrino experiments but one (Homestake) have been calibrated, and their 
experimental responses were found to be in a very good agreement with expectations. The 
argon extraction efficiency of the Homestake detector was also checked by doping it with a 
known small number of radioactive argon atoms, but no calibration has been carried out 
since no artificial source of neutrinos with a suitable energy spectrum exists. 

The solar neutrino problem is not just the problem of the deficit of the observed neutrino 
flux: results of different experiments seem to be inconsistent with each other. In the absence 
of new neutrino physics, the energy spectra of the various components of the solar neutrino 
flux are given by the standard nuclear physics and well known, and only the total fluxes of 
these components may be different from those predicted by the standard solar models. One 
can then infer the flux of neutrinos directly from the Kamiokande or Super-Kamiokande 
data, and use this flux to flnd the corresponding contribution to the Homestake detection 
rate. This contribution turns out to be larger than the total detection rate, so the best-flt 
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Figure 11: Predictions of standard solar models since 1988. The ^Be and fluxes are normalized 
to the predictions of the BP98 model. The rectangular box defines the 3cr error range of the BP98 
model. The best-fit ''Be flux is negative (from ref. p9| ). 
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Figure 12: Solar neutrino measurements and theoretical flux predictions. For CI and Ga ex- 
periments, the units are SNU (1 SNU = 10~^^ captures per target atom per second); for H2O 
experiments, the ratio data/exp. is shown. From ref. ||4^ . 
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contribution of the ^Be neutrinos to the Homestake detection rate is negative! In fact, if one 
assumes that the solar neutrino spectra are undistorted, one can demonstrate the existence 
of the solar neutrino problem without using any information about solar physics except the 
solar luminosity constraint, or even without this constraint, i.e. in a solar-model independent 
way. Therefore the probability of the astrophysical solution of the solar neutrino problem 
is very low. 

Alternatively, one could assume that the observed deficit of solar neutrinos is due to some 
unknown experimental errors. However, in this case it is not sufficient to assume that one 
of the solar neutrino experiments is wrong; one would have to assume that the experiments 
of at least two different types (chlorine, gallium and water Cherenkov) are wrong, which is 
very unlikely. The remaining possibility is that the spectra of solar neutrinos are distorted, 
which requires new neutrino physics. 

There are several possible particle-physics solutions of the solar neutrino problem, the 
most natural one being neutrino oscillations (for an alternative solution related to nonzero 
neutrino magnetic moments, see e.g. [0). The neutrino oscillation solution has become 
even more plausible after the strong evidence for atmospheric neutrino oscillations was re- 
ported by the Super-Kamiokande Collaboration. Neutrino oscillations can convert a fraction 
of solar z/g into u^j, or Ur (or their combination). Since the energy of solar neutrinos is smaller 
than the masses of muons and tauons, these z/^ or Ur cannot be detected in the CC reactions 
of the type ( |129| ) or ( |130|) and therefore are invisible in the chlorine and gallium experiments. 
They can scatter on electrons through the NC interactions and therefore should contribute to 
the detection rate in water Cherenkov detectors. However, the cross section of the NC chan- 
nel of reaction (|131|) is about a factor of 6 smaller than that of the CC channel, and so the 
deficit of the neutrino fiux observed in the Kamiokande and Super-Kamiokande experiments 
can be explained. The probabilities of neutrino oscillations depend on neutrino energy, and 
the distortion of the energy spectra of the experimentally detected solar neutrinos, which is 
necessary to reconcile the data of different experiments, is readily obtained. 

The oscillations of solar neutrinos inside the sun can be strongly enhanced due to the 
MSW effect which we discussed in sec. |], and the solar data can be fitted even with a 
very small vacuum mixing angle. Solar matter can also infiuence neutrino oscillations if the 
vacuum mixing angle is not small. The allowed values of the neutrino oscillation parameters 
sin^ 26q and Am^ which fit the detection rates in the chlorine, gallium and water Cherenkov 
experiments in the 2-fiavour scheme are shown in fig. |T3 for matter enhanced (left panel) 
and vacuum (right panel) oscillations of z/g into active neutrinos. In the case of the matter 
enhanced oscillations, there are three allowed ranges of the parameters corresponding to 
the small mixing angle (SMA), large mixing angle (LMA) and low Am^ (LOW) MSW 



solutions [50, 0. The LOW solution has a very low probability and only exists at 99% 
C.I.. The vacuum oscillation (VO) solution corresponds to very small values of Am^, for 
which the neutrino oscillation length (|6l|) for typical solar neutrino energies (~ a few MeV) 
is comparable to the distance between the sun and the earth. This solution is also known as 
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Figure 13: Left panel: Allowed regions (99 % c.l.) for the MSW solutions of the solar neutrino 
problem for 2-flavour oscillations of i^e into active neutrinos (only rates included). Right panel: 
The same for the VO solution. From ref. pO|]. 

"just so" oscillation solution. For z/g — >■ i^steriie oscillations, there is only the SMA solution 
with the allowed region of parameters similar to that for oscillations into active neutrinos. 

Solar neutrinos detected during night travel some distance inside the earth on their 
way to the detector, and their oscillations can be affected by the matter of the earth. 
In particular, a fraction of or Vj. produced as a result of the solar v^. oscillations can 
be reconverted into Vf. by oscillations enhanced by the matter of the earth. This earth 
"regeneration" (or day/night) effect can be appreciable in the case of the LMA solution, 
but is expected to be very small in the case of the SMA solution. The day/night effects (and 
in general, the zenith angle dependence of the neutrino signal) can in principle be observed 
in real-time experiments, such as Super-Kamiokande. 

In the case of the MSW solutions of the solar neutrino problem the neutrino state 
arriving at the earth is an incoherent superposition of the mass eigenstates vi and V2 (see 
Problem 23 below). The probability Pse of finding a solar v^. after it traverses the earth 
depends on the average Ve survival probability in the sun Pg: 

1 op 

Pse = Ps + ^ {P2e - sin^ ^o) • (132) 

cos 26'o 

Here = P{^2 ~^ ^e) is the probability of oscillations of the second mass eigenstate into 
electron neutrino inside the earth. It follows from this expression that the sign of the night- 
day asymmetry {N ~ D) / {N -\- D) depends on whether Ps is larger or smaller than 1/2. For 
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the LMA solution, Ps < 1/2, and the night-day asymmetry is expected to be positive. For 
the SMA solution, Ps for the Super-Kamiokande experiment is close to 1/2, and therefore 
the expected day/night effect is small. It can in principle be of either sign. The value 
measured in the Super-Kamiokande experiment is ||5^ 



A^ 

_ - 1 = 0.067 ± 0.033 (stat.) ± 0.013 (syst.) , (133) 

i.e. shows an excess of the night-time flux at about 2a. The zenith angle event dependence 
measured by the Super-Kamiokande shows a rather fiat distribution of the excess of events 
over different night-time zenith angle bins. This is rather typical for the LMA solution 
whereas for the SMA solution one normally expects the excess (or deficiency) to be con- 
centrated in the vertical upward bin with zenith angles 6 in the range —1 < cosO < —0.8. 
Thus, the night-time zenith angle dependence and the overall night-day asymmetry ( |133| ) 
favour the LMA solution of the solar neutrino problem. 

As we discussed before, oscillations of solar neutrinos in the 3-fiavour scheme can be 
described through the effective two-flavour z/g survival probability (see eq. (|7D and Prob- 
lem 20). For a recent detailed analysis of the solar neutrino oscillations in the 3-fiavour 
framework, see [p3 |. 



Problem 22. Discuss effects of the solar and earth 's matter on the oscillations of solar 
neutrinos in the case of VO solution. Take into account that the matter densities psun o,nd 
pearth Vary in the ranges - 150 and 2.7 - 12.5 g/cm? respectively. 

Problem 23. Show that in the case of the MSW solutions the solar neutrino flux arriving 
at the earth can be considered as an incoherent superposition of the neutrino mass eigenstates 
Vi and i>2- Hint: consider the relative phase between vi and V2 states acquired by neutrinos 
on their way from the sun to the earth. Assume that the energy resolution of the neutrino 
detector AE > 100 keV. 



Problem 24. Derive equation i\132l) for the probability of finding a solar after it crosses 



the earth during night. Hint: use the results of the previous problem. 

As we have already discussed, neutrino oscillations (both in matter and in vacuum) 
depend on neutrino energy and therefore should result in certain distortions of the spectrum 
of the detected solar neutrinos. One can obtain an information on the spectrum of incoming 
neutrinos by measuring the recoil electron spectrum in the reaction (|131|) . This spectrum has 
been measured in the Super-Kamiokande experiment, and the results (along with the typical 
predictions of the LMA, SMA and VO solutions) are shown in fig. ([I^ . In the absence of the 
neutrino spectrum distortion, the ratio of measured/expected electron spectrum presented in 
this figure should be a horizontal line. The characteristic feature of the measured spectrum is 
an excess of events in the region of high energies. The experimental errors are still too large 
to allow a clean discrimination between different solutions of the solar neutrino problem; 
however, the spectrum favours VO over the other solutions. Another possible explanation 
of the excess of the high energy events is that the flux of the highest-energy hep neutrinos 
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Figure 14: Recoil electron energy spectrum in the Super-Kamiokande experiment along with 
typical predictions of SMA, LMA and VO solutions (from ref. 



has been underestimated by about a factor of 20. Since the hep flux is very poorly known, 
this is still a possibility. 

At present, the situation with the solar neutrino data is somewhat confusing: different 
pieces of the data favour different solutions of the solar neutrino problem. The total rates 
are better fitted within the SMA solution, the night/day ratio prefers the LMA solution, 
and the electron recoil spectrum in Super-Kamiokande favours the VO solution. Clearly, 
more data are needed to clear the situation up. Fortunately, two new experiments which 
can potentially resolve the problem are now under way or will soon be put into operation. 
The SNO (Sudbury Neutrino Observatory) experiment has started taking data, and its first 
results are expected some time during 2000. The Borexino experiment is scheduled to start 
data taking in 2001. 

The SNO detector consists of 1000 tons of heavy water, and it detects solar neutrinos 
in three different channels: 

u^ + d^p + p + e- (CC), E™„ = 1.44MeV, (134) 
Ua + d^p + n + ua (NC), E™„ = 2.23MeV, (135) 

and scattering process ( |131| ) which can proceed through both CC and NC channels. The 
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CC reaction (|134|) is very well suited for measuring the solar neutrino spectrum. Unlike in 



the case of scattering ( |131| ) in which the energy of incoming neutrino is shared between 
two light particles in the final state, the final state of the reaction ( |134| ) contains only 
one light particle - electron, and a heavy 2p system whose kinetic energy is very small. 
Therefore by measuring the electron energy one can directly measure the spectrum of the 
solar neutrinos. Detecting the solar neutrino spectrum would be of great importance since 
any deviation from the spectrum predicted by the nuclear beta decay would be a "smoking 
gun" signature of new neutrino physics. 

The cross section of the NC reaction (|135|) is the same for neutrinos of all three flavours. 



and therefore oscillations between i/^ and or Ur would not change the NC detection rate in 
the SNO experiment. On the other hand, these oscillations would deplete the solar z/g flux, 
reducing the CC event rate. Therefore the CC/NC ratio is a sensitive probe of neutrino 
flavour oscillations. If solar Uf, oscillate into electroweak singlet (sterile) neutrinos, both CC 
and NC event rates will be suppressed. 

The Borexino experiment will detect solar neutrinos through the scattering with 
a very low energy threshold, and will be able to detect the ^Be neutrino line. Different 
solutions of the solar neutrino problem predict different degree of suppression of ''Be neu- 
trinos, and their detection could help discriminate between these solutions. Observation of 
the ^Be neutrino line would be especially important in the case of the VO solution. Due to 
the eccentricity of the earth's orbit the distance between the sun and the earth varies by 
about 3.5% during the year, and this should lead to varying oscillation phase (and therefore 
varying solar neutrino signal) in the case of vacuum neutrino oscillations. This seasonal 
variation can in principle be separated from the trivial 7% variation due to the law 
which is not related to neutrino oscillations. However, the oscillation phase depends on 
neutrino energy [see (^9|)], and integration over significant energy intervals may make it 
difficult to observe the seasonal variations of the solar neutrino flux due to the VO. The 
^Be neutrinos are monochromatic, which should facilitate the observation of the seasonal 
variations. The oscillation phase in ( |5DD depends on the ratio L/E, therefore for neutrinos 
with continuous spectrum VO should lead to seasonal variations of the spectrum distortion 



which are correlated with the seasonal variations of the flux |Q. Such variations would be 
a clear signature of the VO solution. 

Although up to now the sensitivities of the solar neutrino experiments have been insuf- 
ficient to allow a clear-cut discrimination between different solutions of the solar neutrino 
problem, one can hope that the combined data of the currently operating and forthcoming 
experiments will allow to finally resolve the solar neutrino problem. 
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11 Atmospheric neutrinos 



Atmospheric neutrinos are electron and muon neutrinos and their antineutrinos which are 
produced in the hadronic showers induced by primary cosmic rays in the earth's atmosphere. 
The main mechanism of production of the atmospheric neutrinos is given by the following 
chain of reactions: 

p{a,...) + Air ir^iK^) + X 



IT 



±(K±) ^ fi^ + u^iu^) (136) 



Atmospheric neutrinos can be observed directly in large mass underground detectors pre- 
dominantly by means of their CC interactions: 

Pe{ye) + A^ e-(e+) + X, 

z/^(z7^) + A^^-(^+)+X. (137) 



Naively, from the reaction chain (|136|) one would expect to have two atmospheric muon neu- 
trinos or antineutrinos for every electron neutrino or antineutrino. In reality, the situation 
is more complicated: one should take into account the differences in the lifetimes of vr^, 
and /i^ as well as the differences in their spectra. Also, although the reaction chain ( |136| ) 
is the dominant source of atmospheric neutrinos, it is not the only one. Calculations of the 
atmospheric neutrino fluxes predict the v^^/ Ve ratio that depends on neutrino energy and the 
zenith angle of neutrino trajectory, approaching 2 for low energy neutrinos and horizontal 
trajectories but exceeding this value for higher energy neutrinos and for trajectories close 
to vertical. 

Accurate calculation of the atmospheric neutrino fluxes is a difficult job which includes 
such ingredients as spectra and chemical composition of cosmic rays (including geomag- 
netic effects and solar activity), cross sections of tt and K production off the nuclear tar- 
gets, Monte Carlo simulation of hadronic cascades in the atmosphere and the calculation 
of neutrino spectra including muon polarization effects. Each step introduces some un- 
certainty in the calculation. The overall uncertainty of the calculated atmospheric neu- 
trino fluxes is rather large, and the total fluxes calculated by different authors differ by as 
much as 20 - 30%. At the same time, the ratio of the muon to electron neutrino fluxes 
is fairly insensitive to the above uncertainties, and different calculations yield the ratios 
of muon-like to electron-like contained events which agree to about 5%. This ratio has 
been measured in a number of experiments, and the Kamiokande and 1MB Collabora- 
tions reported smaller than expected ratio in their contained events, with the double ratio 
R{fi/e) = [(z/^ + + ^^e]data/[{l^^l + ^ti)/{.^e + i^eTlMC - 0.6 where MC stands for Monte 

Carlo simulations. The discrepancy between the observed and predicted atmospheric neu- 
trino fluxes was called the atmospheric neutrino anomaly. The existence of this anomaly 
was subsequently confirmed by Soudan 2, MACRO and Super-Kamiokande experiments. 
Most remarkably, the Super-Kamiokande (SK) Collaboration obtained a very convincing 
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evidence for the up-down asymmetry and zenith-angle dependent deficiency of the flux of 
muon neutrinos, which has been interpreted as an evidence for neutrino oscillations. We 
shall now discuss the SK data and their interpretation. 

The SK detector is a 50 kt water Cherenkov detector (22.5 kt fiducial volume) which 
is overseen by more than 13,000 photomultiplier tubes. The charged leptons born in the 
CC interactions of neutrinos produce the rings of the Cherenkov light in the detector which 
are observed by the phototubes. Muons can be distinguished from electrons since their 
Cherenkov rings are sharp whereas those produced by electrons are diffuse. The SK Collab- 
oration subdivided their atmospheric neutrino events into several groups, depending on the 
energy of the charged leptons produced. Fully contained (FC) events are those for which 
the neutrino interaction vertex is located inside the detector and all final state particles do 
not get out of it. FC events are further subdivided into sub-GeV (visible energy < 1.33 
GeV) and multi-GeV (visible energy > 1.33 GeV) events. Partially contained (PC) events 
are those for which the produced muon exits the inner detector volume (only muons are 
penetrating enough). The average energy of a neutrino producing a PC event in SK is ~ 15 
GeV. Muon neutrinos can also be detected indirectly by observing the muons that they 
have produced in the material surrounding the detector. To reduce the background from 
atmospheric muons, only upward-going neutrino-induced muons are usually considered. A 
rough estimate of the energy spectrum of the upward-going muons has been obtained di- 
viding them in two categories, passing (or through-going) and stopping muons. The latter, 
which stop inside the detector, correspond to the average parent neutrino energy ~ 10 GeV, 
whereas for the through-going muons the average neutrino energy is ~ 100 GeV. 

The measurements of the double ratio R{fj,/e) for contained events at SK (848 hve days) 
give H 

R = 0.68 ± 0.02 (stat.) ± 0.05 (syst.) (sub-GeV) , 

R = 0.68 ±OM{stat.)± 0.08 (syst.) (multi-GeV) . (138) 

The value of R for sub-GeV events is different from unity (to which it should be equal in 
no-oscillation case) by 5.9a. 

We shall now discuss the zenith angle distributions of the atmospheric neutrino events. 
It should be remembered that the zenith angle distributions of the charged leptons which 
are experimentally measured do not coincide with those of their parent neutrinos: for multi- 
GeV neutrinos the average angle between the momenta of neutrinos and charged leptons 
is about 17°, whereas for sub-GeV neutrinos it is close to 60°. This is properly taken into 
account in MC simulations. For PC events and upward going muons the correlation between 
the directions of momenta of muons and parent neutrinos is much better. The distances 
L traveled by neutrinos before they reach the detector vary in a wide range: for vertically 
downward going neutrinos (neutrino zenith angle = 0) L ~ 15 km; for horizontal neutrino 
trajectories {6^, = 90°) L ~ 500 km; the vertically up-going neutrinos {9u = 180°) cross the 
earth along its diameter and for them L ~ 13, 000 km. 
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Figure 15: Zenith angle distributions for sub-GeV and multi-GeV e-like and ^-like events at SK. 
The bars show the (no-oscillations) Monte Carlo predictions; the lines show the predictions for 
i^fi ^ i^T oscillations with the best-fit parameters Am^ = 3.5 x 10~^ eV^, sin^ 29 = 1.0. From [^5| . 



In fig. 15 the zenith angle distributions of the SK e-like and /i-like events are shown 



separately for sub-GeV and multi-GeV contained events. The data correspond to 736 live 
days. (By the time these lecture notes were being written, the results for 848 live days of 
data taking had been reported, but no high-quality figures for zenith angle distributions 
were publicly available). One can see that for e-like events, the measured zenith angle dis- 
tributions agree very well with the MC predictions (shown by bars), both in the sub-GeV 
and multi-GeV samples, while for /i-like events both samples show zenith-angle dependent 
deficiency of event numbers compared to expectations. The deficit of muon neutrinos is 
stronger for upward going neutrinos which have larger pathlengths. In the multi-GeV sam- 
ple, there is practically no deficit of events caused by muon neutrinos coming from the upper 
hemisphere (cos 6' > 0), whereas in the sub-GeV sample, all /i-like events exhibit a deficit 
which decreases with cos^. This pattern is perfectly consistent with oscillations z/^ ^ u^- 
or z/^ i>s where is a sterile neutrino. Muon neutrinos responsible for the multi-GeV 
sample are depleted by the oscillations when their pathlength is large enough; the deple- 
tion becomes less pronounced as the pathlength decreases (cos 6' increases); for neutrinos 
coming from the upper hemisphere, the pathlengths are too short and there are practically 
no oscillations. Neutrinos responsible for the sub-GeV /x-like events have smaller energies. 
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and so their oscillation lengths are smaller; therefore even neutrinos coming from the up- 
per hemisphere experience sizeable depletion due to the oscillations. For up-going sub-GeV 
neutrinos the oscillation length is much smaller than the pathlength and they experience 
averaged oscillations. The solid line in fig. |T5| obtained with the v-r oscillation param- 

eters in the 2-flavour scheme Am^ = 3.5 x 10~^ eV^, s\.v?26 = 1.0 gives an excellent fit of 
the data. 
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Figure 16: Left panel: up-down asymmetry vs event momentum for single ring e-like 
events in SK. Right panel: SK allowed regions for oscillations parameters for 
2-flavour scheme (FC -|- PC single ring events). From Q. 
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An informative parameter characterizing the distortion of the zenith angle distribution 
is the up-down asymmetry A, where up corresponds to the events with cos^ < —0.2 and 
down to those with cos^ > 0.2. The flux of atmospheric neutrinos is expected to be nearly 
up-down symmetric for neutrino energies E > 1 GeV, with minor deviations coming from 
geomagnetic effects which are well understood and can be accurately taken into account. 
In particular, at the geographical location of the SK detector small positive asymmetry is 
expected. Any significant deviation of the up-down asymmetry of neutrino induced events 
from the asymmetry due to the geomagnetic effects is an indication of neutrino oscillations 
or some other new neutrino physics. The asymmetry measured for the SK multi-GeV /^-like 
events is ^ 

A = = -0.32 ± 0.04 (stat.) ± 0.01 (syst.) , (139) 

i.e. is negative and differs from zero by almost Sal The dependence of the asymmetries for 



e-like FC and /i-like FC-I-PC events on the event momentum is shown in fig. |T6| (left panel). 
One can see that for e-like events A ~ for all momenta. At the same time, for /i-like events 
the asymmetry is close to zero at low momenta and decreases with momentum. This is easily 
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understood in terms of the i/^ oscillations. For very small momenta, the oscillation length 
is small and both up-going and down-going neutrino fluxes are depleted by oscillations to 
about the same extent; in addition, loose correlation between the directions of the momenta 
of the charged lepton and of its patent neutrino tends to smear out the asymmetry at low 
energies. With increasing momentum the oscillation length increases, and the pathlength 
of down-going neutrinos becomes too small for oscillations to develop. The right panel of 
fig. O gives the allowed values of the oscillation parameters following from the SK FC and 



PC event data. The best fit corresponds to Am = 3.1 x 10~ eV , sin 29 = 0.99 and has a 
very good /d-o.f. = 55/67. In this analysis the poorly known overall normalization factor 
of the neutrino flux was considered as a free parameter, and the best fit was achieved with 
5% upward renormalization of the flux of Honda et al. |5^. The above value of should 
be compared with that of the no-oscillation hypothesis: /d.o.f. = 177/69, which is a very 
poor fit. 
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Figure 17: Left panel: SK allowed regions of oscillation parameters for z/^ channel in 2-flavour 
scheme (FC + PC single ring events + upward through-going -|- upward stopping events). Right 
panel: comparison of allowed values of oscillation parameters obtained in different experiments. 
From ||. 

The SK data show evidence for neutrino oscillations not only in their FC and PC /x-like 
events: upward stopping and upward through-going events also demonstrate zenith angle 
dependent deficiency of muon neutrinos consistent with neutrino oscillations, although the 
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statistics for up-going muons is lower than that for contained events. Fig. (left panel) 
summarizes the allowed region for 2-flavour neutrino oscillations obtained from the combined 
analysis of the SK FC and PC single ring events (848 live days), upward through-going (923 
live days) and upward stopping events (902 live days). The best fit {x^/d.o.f = 67.5/82) was 
obtained for Am^ = 3.5 x 10^^ eV^, sin^ 29 = 1.0 and the upward 6% renormalization of the 
flux of Honda et al. The no-oscillation fit gives x^/d-o.f. = 214/84. The right panel gives the 
comparison of the allowed regions obtained in different atmospheric neutrino experiments. 
The evidence for neutrino oscillations is very impressive. 

Are neutrino oscillations that are responsible for the depletion of the z/^ flux z/^ Uj. 
or Ufj, ^ Ugl For contained events, the oscillation probabilities in these two channels are 
nearly the same and the data can be fitted equally well in both cases, with very similar 
allowed ranges of the oscillation parameters. However, for higher energy upward going events 
there are important differences between these two cases. In the 2-flavour scheme, z/^ ^ 
oscillations are not affected by matter because the interactions of z/^ and z/^ with matter 
are identical. However, sterile neutrinos do not interact with matter at all, and therefore 
the <-> Us oscillations are affected by the matter- induced potential — Vg = V^. At low 
energies, the kinetic energy difference Am? /2E dominates over V^, and the earth's matter 
effects are unimportant. They become important at higher energies, when An? /2E ~ 



at very high energies, when Am /2E -C V^, matter strongly suppresses neutrino oscillations 
both in <-> Vg and ^ channels (see eq. (|97| ) where has to be replaced by =fA^„/2 
for these oscillation channels). Therefore the oscillations of high energy neutrinos travelling 
significant distances in the earth should be strongly suppressed in this case. 

Fig. |l^ (upper left panel) shows the zenith angle distribution of the upward through- 
going events {{E) ~ 100 GeV) along with the predictions of z/^ ^ Ur and z/^ ^ z/^ scenarios 
with parameters obtained from the fits of contained events. It can be seen that the z/^ ^ Ur 
oscillation channel describes the "vertical" (cos^ < —0.4) events better than the ^ Ug 
channel. For horizontal and nearly horizontal events (—0.4 < cos 6' < 0) both scenarios give 
similar predictions since the neutrinos do not travel inside the earth or travel only small 
distances there. This is further illustrated by the left bottom panel of fig. |18| in which the 
Vertical/Horizontal event ratio is plotted. The excluded ranges of the oscillation parameters 



for various oscillation channels are shown in the right panel of fig. |T8|. One can conclude 
that the z/^ Ug oscillations are currently disfavoured at the 2a level, although not yet 
completely ruled out as a possible solution of the atmospheric neutrino problem. 

Another way to discriminate between z/^ Ur and z/^ ^ Ug as the two possible main 
channels of the atmospheric neutrino oscillations is through the detection of single pions 
produced in the NC reactions uN —>■ uNn^. Sterile neutrinos do not participate in these 
reactions and therefore Ug oscillations would deplete the number of the observed 7r*^'s. 

Unfortunately, cross sections of these NC reactions are only known with large uncertainties, 
and at the moment a clear discrimination between the z/^ and z/^ Ug oscillations 

of atmospheric neutrinos through the detection of single neutral pions does not seem to be 
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Figure 18: Left top panel: Zenith angle distribution of upward through-going events at SK. The 
dashed and short-dashed histograms show the Vr and i/^ i/g solutions respectively. Left 

bottom: Vertical/Horizontal ratio for upward through-going muons along with the predictions of 
i/^ <-> Ur and ^ Us oscillation scenarios. Right panel: Exclusion regions obtained from Up/Down 
ratio for energetic PC sample {Ems > 5 GeV, {E) ~ 25 GeV) and Vertical/Horizontal ratio for 
upward through-going events. From |^]. 



possible. The situation may change in the near future, however, as the K2K Collaboration 
is planning to perform high statistics measurements of neutrino interactions in their near 
detector. 

Can ^ Vf. oscillations be responsible for the observed anomalies in the atmospheric 
neutrino data? The answer is no, at least not as the dominant channel. Explaining the data 
requires oscillations with large mixing equal or close to the maximal one; v^. oscillations 

would then certainly lead to a significant distortion of the zenith angle distributions of the 
e-like contained events, contrary to observations. In addition, for Am^ in the range ~ 10"'^ 
eV^ which is required by the atmospheric neutrino data, z/^ ^ Vf, oscillations are severely 
restricted by the CHOOZ reactor antineutrino experiment, which excludes these oscillations 
as the main channel of the atmospheric neutrino oscillations (see sec. |12D. 



However, z/„ ^ z/g and z/g ^ Vt can be present as subdominant channels of the oscilla- 
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tions of atmospheric neutrinos. This may lead to interesting matter effects on oscillations 
of neutrinos crossing the earth on their way to the detector. Matter can strongly affect 
J^e ^ i^ti,T and Ue ^ i^fi,T oscillations leading to an enhancement of the oscillation probabili- 
ties for neutrinos and suppression for antineutrinos or vice versa, depending on the sign of 
the corresponding mass squared difference. For certain neutrino energies, one can expect an 
enhancement due to the MSW effect in the mantle or in the core of the earth. In addition, 
there can be a specific enhancement effect, different from the MSW one, due to a cumu- 
lative build-up of the transition probability for neutrinos crossing the earth's mantle, core 



and then again mantle (parametric resonance ||4^, ^ 0)- The matter enhancement 
effects are expected to take place in a relatively narrow range of neutrino energies, and it 
is rather difficult to observe them experimentally. Yet the possibility of the observation of 
matter effects on oscillations of atmospheric neutrinos is very exciting. 

Since three neutrino flavours are known to exist, oscillations of atmospheric neutrinos 
should in general be considered in the full 3-flavour framework (assuming that no sterile 
neutrinos take part in these oscillations). As follows from the analyses of solar and at- 
mospheric neutrino data, there should be two distinct mass squared difference scales in 
the three neutrino framework, Am^^^ ~ lO^'^eV^ and Am|, < 10~^eV^. The hierarchy 
Am^^^ ^ ArriQ means that one of the three neutrino mass eigenstates (which we denote 
1/3) is separated by the larger mass gap from the other two. The neutrino mass squared dif- 
ferences thus satisfy |Am3ip ~ |Am32p ^ |Am2ip. We know already that z/^ should 
be the main channel of oscillations whereas ^ v^^t oscillations can only be present as 
the subdominant channels. In the 3-fiavour framework this means that the element \Uez\ 
of the lepton mixing matrix (0) is small. The smallness of \Uez\ means that the <-> Vj. 
and Pe oscillations approximately decouple, and the 2-fiavour analysis gives a good 
first approximation. In terms of the parametrization ( p^D of the lepton mixing matrix, the 
mixing angle describing the main channel of the atmospheric neutrino oscillations is ^23, 
and its best-fit value following from the SK data is 45°. This fact and the smallness of 
If^esl = 513 mean that the mass eigenstate 1/3 mainly consists of the flavour eigenstates 
and Vj. with approximately equal weights, while the admixture of to this state is small or 
zero. 

Can one observe true 3-fiavour effects in the atmospheric neutrino oscillations? In 
principle, yes. Although in the 2-fiavour scheme v-r oscillations are not affected by 

matter, in the full 3-fiavour approach these oscillations are modified by matter because of 
the mixing with z/g. These matter effects are, however, rather small because of the smallness 
of \Uez\- Three- flavour effects can also lead to interesting phenomena in the subdominant 

^ and z/g ^ oscillations. To see this, let us consider the fluxes and of 



atmospheric z/^ and z/g in the presence of oscillations. As we discussed in sec. |8.4| , in the 
case of the hierarchy |Am3ip ~ |Am32p ^ |Am2ip the oscillation probabilities in the 3- 
flavour scheme simplify significantly and can be expressed through the evolution matrix of 
the effective 2-flavour problem. Using eqs. ( |1 19| ) - ( |123| ) one can express F^ and Fg through 
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the corresponding no-oscillation fluxes F° and as 



F, = F0[1 + P2M3-1)] > (140) 



F = F^ 



1 - % 1 + 2sl^4^ (Zcos<l> + l^asin^ - 1] 



where r{E,9y) = Fl^{E,9^)/F^{E,9^). It is interesting that the subdominant 
oscillations may lead either to an excess or to a deficiency of e-like events depending on 
the values of r and S23. Indeed, the effect of oscillations on the electron neutrino flux is 
proportional to the factor (rs23 ~ 1)- 1^ assumes r = 2, there will be an excess of 
e-like events for 6*23 > 45° and a deficiency for 6*23 < 45°. The SK best fit is 6*23 = 45°; 
for this value of ^23 there would be no deviation from the prediction for r = 2. However, 
for upward going neutrinos in the multi-GeV range r is typically 3 - 3.5 rather than 2, so 
there should be an excess of e-like events for 6*23 > 33° and a deficiency for 6*23 < 33°. Thus, 
the distortions of the zenith angle distributions of the e-like events due to the subdominant 
^ ^^l,T oscillations should depend on the value of the mixing angle 6*23 governing the 
dominant ^ oscillation channel. However, these distortions are expected to be small 
because of the smallness of |[/e3| (notice that the probability in ( |140| ) and ( |141[ ) vanishes 
in the limit f/e3 0). For recent analyses of the SK atmospheric neutrino data in the 
3-flavour framework, see IF? 



If the solution of the solar neutrino problem is the LMA MSW effect, lS.vn\^ can be as 
large as ~ lO""^ eV^ and its influence of the atmospheric neutrino oscillations cannot be 
neglected. In that case eqs. (|140D and (|141D have to be modified. It was shown in |SH] that 



in this case the subdominant oscillations can lead to a negative up-down asymmetry of e- 
like events at relatively large momenta. If the trend of slightly decreasing with momentum 
asymmetry of the e-like events which can be seen in the left top panel of fig. |16| is of 
physical origin and not just due to statistical fluctuations, it can be explained by relatively 
large lS.v!i^^. 

Are the standard neutrino oscillations the sole possible explanation of the observed at- 
mospheric neutrino anomalies? In principle, other explanations are possible. Those include 
exotic types of neutrino oscillations - matter-induced oscillations due to flavour- changing 



interactions of neutrinos with medium |£9[, oscillations due to small violations of the Lorentz 
or CPT invariance or of the gravitational equivalence principle , and also neutrino decay 
6T| . Exotic oscillations lead to periodic variations of the v^^ survival probability with the 



oscillation lengths lose oc where n = in the case of flavour-changing neutrino interac- 
tions or violation of CPT invariance, and n = 1 for oscillations due to the violations of the 
Lorentz invariance or equivalence principle. This has to be contrasted with n = — 1 in the 
case of the standard neutrino oscillations. The energy dependence of the oscillation length 
can be tested in the atmospheric neutrino experiments as the energies of detected neutrinos 



span more than 3 orders of magnitude. The analysis was performed in and the authors 
found that the fit of the atmospheric neutrino data assuming oscillations with Iq^c oc 
gives n = —0.9 ± 0.4 at 90% c.l.. This clearly favours the standard oscillations over the 
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exotic ones. In contrast to this, the neutrino decay mechanism fits the SK data quite well, 
the quality of the fit being as good as the one for the standard neutrino oscillations [ p^ . 
However, the assumed neutrino decay mode (into a sterile neutrino and a Majoron) requires 
particle physics models which look less appealing than those that provide neutrino masses 
required for the standard neutrino oscillations. 

Problem 25. Using eqs. ^lldj) - ( \12!^ ) of sec. \8.4\ derive expressions ( \14(] ) and ( \141\ ) for 
the atmospheric z/g and fluxes in the presence of oscillations in the case of the hierarchy 
|Am3ip| ~ Am32p ^ |Am2ip. 



12 Reactor and accelerator neutrino experiments 

In reactor neutrino experiments oscillations of electron antineutrinos into another neutrino 
species are searched for by studying possible depletion of the z/g flux beyond the usual 
geometrical one These are the disappearance experiments, because the energies 

of the reactor z/g's {{E) ~ 3 MeV) are too small to allow the detection of muon or tauon 
antineutrinos in CC experiments. Small energy makes the reactor neutrino experiments 
sensitive to oscillations with rather small values of Arn^. 
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Figure 19: Left panel: CHOOZ exclusion region for z^e z^x oscillations. The region to the 
right of the curves are excluded. Shaded areas correspond to the Kamiokande allowed region for 
atmospheric z^^ — > Ve oscillations |6(:]. Right panel: LSND allowed region for — > z^g oscillations 
(shaded areas) along with KARMEN, BNL776 and Bugey exclusion regions 



Up to now, no evidence for neutrino oscillations has been found in the reactor neutrino 
experiments, which allowed to exclude certain regions in the neutrino parameter space. 
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The best constraints were obtained by the CHOOZ experiment in France |^ (see fig. 
|TU|). As follows from eq. ([7^), in the case of the hierarchy Amg^^ ~ Am^g ^ Am2i, 
the survival probability of Ve (which by CPT invariance is equal to that of Ue) in the 3- 
fiavour scheme has a simple form which coincides with that in the two-fiavour framework. 
Therefore the excluded regions obtained in the 2-fiavour analysis apply to the 3-fiavour 
case as well, with the identification Am^ = Am|^ ~ Am^j, sin 26 = sin 26'i3. For the 
values of Amg^ = Am^^^ in the Super-Kamiokande allowed region (2 — 6) ■ 10~^ eV^, 
the CHOOZ results give the following constraint on the element Ues of the lepton mixing 
matrix: |[/e3p(l — |t/e3p) < 0.045 — 0.02, i.e. \Ue3\ is either small or close to unity. The 
latter possibility is excluded by solar and atmospheric neutrino observations, and one finally 
obtains 

sin^ ^13 = |t/e3p < (0.047 - 0.02) for Am^^ = (2 - 6) ■ 10"^ eV^ . (142) 
This is the most stringent constraint on {Uesl date. 

Presently, a long baseline reactor experiment KamLAND is under construction in Japan. 
This will be a large liquid scintillator detector experiment using the former Kamiokande site. 
KamLAND will detect electron antineutrinos coming from several Japanese power reactors 
at an average distance of about 180 km. KamLAND is scheduled to start taking data in 
2001 and will be sensitive to values of Am^ as low as 4 x 10~^ eV^, i.e. in the range relevant 
for the solar neutrino oscillations! It is expected to be able to probe the LMA solution of 
the solar neutrino problem (see fig. ^0|). It may also be able to directly detect solar 
and ''Be neutrinos after its liquid scintillator has been purified to ultra high purity level by 
recirculating through purification p7[ |. 



In the 3-fiavour framework, the oscillations due to the large Am^^^ ~ 3 x 10~^ eV^ 
should be in the averaging regime in the KamLAND experiment and therefore the survival 
probability should be described by eqs. (|74|) and (|75| ) rather than by ([72|). In the case of 



SMA, LMA and VO solutions of the solar neutrino problem the z/g survival probability at 
KamLAND is then just P{ue ^e) — C13, i.e. the mixing angle 613 is directly measured. The 
CHOOZ data allow > 0.90 for Am^^^ in the SK range. In the case of the LMA solution. 
Pipe ^e)zfi — c.\.^Pipe ^e)2fi where Pipe ^e)2fi IS the usual 2-fiavour survival 
probability (|75D. Thus, in this case the 2-fiavour analysis is modified by a non- vanishing 
On- 

There have been a number of accelerator experiments looking for neutrino oscillations. 
In all but one no evidence for oscillations was found and constraints on oscillation parameters 



were obtained [|65| , p8|] . The LSND Collaboration have obtained an evidence for z/^ — > z/g 



and z/^ — > z/g oscillations [Q. The LSND result is the only indication for neutrino oscillations 
that is a signal and not a deficiency. The KARMEN experiment is looking for neutrino 
oscillations in z/^ — > Vf. channel. No evidence for oscillations has been obtained, and part 
of the LSND allowed region has been excluded. In fig. |19| (right panel) the results from 
LSND and KARMEN experiments are shown along with the relevant constraints from the 
BNL776 and Bugey experiments. One can see that the only domain of the LSND allowed 
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Figure 20: Left panel: results of the present and sensitivities of future fe <-> oscillations searches 
(90% c.l.) [p^. Right panel: sensitivities of MINOS and K2K long baseline experiments p^]. 



region which is presently not excluded is a narrow strip with sin^ 29 ~ 1.5 x 10 ^ — 4x10 ^ 
and Am2 ~ 0.2 - 2 eV^. 

The existing neutrino anomalies (solar neutrino problem, atmospheric neutrino anomaly 
and the LSND result), if all interpreted in terms of neutrino oscillations, require three 
different scales of mass squared differences: AuIq < 10~^ eV^, Am^^^ ~ 10~^ eV^ and 
AmLSND ~ 0-2 eV^. This is only possible with four (or more) light neutrino species. The 
fourth light neutrino cannot be just the 4th generation neutrino similar to u^, and 
because this would be in conflict with the experimentally measured width of boson [see 
eq. (0)]. It can only be an electroweak singlet (sterile) neutrino. Therefore the LSND result, 
if correct, would imply the existence of a light sterile neutrino. 

Out of all experimental evidences for neutrino oscillations, the LSND result is the only 
one that has not yet been confirmed by other experiments. It is therefore very important to 
have it independently checked. This will be done by the MiniBooNE (first phase of BooNE) 
experiment at Fermilab. MiniBooNE will be capable of observing both v^. appearance 
and z/^ disappearance. If the LSND signal is due to — > v^. oscillations, MiniBooNE is 
expected to detect an excess of ~1500 z/g events during its first year of operation, establishing 
the oscillation signal at ~ 8cr level (see fig. EDI). If this happens, the second detector will be 



installed, with the goal to accurately measure the oscillation parameters. MiniBooNE will 
begin taking data in 2002. 
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A number of long baseline accelerator neutrino experiments have been proposed to date. 
They are designed to independently test the oscillation interpretation of the results of the 
atmospheric neutrino experiments, accurately measure the oscillation parameters and to 
(possibly) identify the oscillation channel. The first of these experiments, K2K (KEK to 
Super-Kamiokande) already started taking data in 1999. It has a baseline of 250 km and 
is looking for z/^ disappearance. K2K should be able to test practically the whole region 
of oscillation parameters allowed by the SK atmospheric neutrino data except perhaps the 
lowest- Am^ part of it (see fig. ^). Two other long baseline projects, NuMI - MINOS 
(Fermilab to Soudan mine in the US) and CNGS (CERN to Gran Sasso in Europe), each 
with the baseline of about 730 km, will be sensitive to smaller values of Am^ and should be 
able to test the whole allowed region of SK (fig. |2y). MINOS will look for z/^ disappearance 
and spectrum distortions due to oscillations. It may run in three different energy 

regimes - high, medium and low energy. MINOS is scheduled to start data taking in the 
end of 2002. CERN to Gran Sasso will be an appearance experiment looking specifically for 
^'^t ^ oscillations. It will also probe z/g appearance. At the moment, two detectors 

have been approved for participation in the experiment - OPERA and ICANOE. The whole 
project was approved in December of 1999 and the data taking is planned to begin on May 
15, 2005 §. 

Among widely discussed now future projects are neutrino factories - muon storage rings 
producing intense beams of high energy neutrinos. In addition to high statistics studies 
of neutrino interactions, experiments at neutrino factories should be capable of measuring 
neutrino oscillation parameters with high precision and probing the subdominant neutrino 
oscillation channels, matter effects and CP violation effects in neutrino oscillations 



13 Phenomenological neutrino mass matrices 

The information on neutrino masses and lepton mixing obtained in solar and atmospheric 
the neutrino experiments can be summarized as follows: 

SMA: Am| ~ (4 - 10) ■ 10-*^ eV^ sin^ 2^© ~ (0.1 - 1.0) ■ 10^^ 



LMA 
VO 
Atm 



Am| ~ (2 - 20) ■ 10-5 eV^ , sin^ 29q ~ 0.65 - 0.97 



^0 

Am| ~ (0.5 - 5) ■ 10-1° eV^ , sin^ 26^ ^ 0.6 - 1.0 

Am2^„ ~ (2 - 6) ■ 10-3 eV^ , sin^ 2^,^^ ~ 0.82 - 1.0 (143) 



In the 3-fiavour framework (neglecting the LSND result). Am© has to be identified with 
Am2i, Amatm with /S.m\i ~ Am^j, Oq with 6'i2, and 9 atm with ^23- The 2-flavour analysis 
gives a good first approximation to the results of the 3-flavour studies ^ because the 
mixing angle ^13 is constrained to be rather small by the CHOOZ data ( |142| ). 

A very concise and illuminating way of summarizing the neutrino data is to present 
the neutrino mass eigenstates graphically as rectangles of unit area, the position of which 
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Figure 21: 3-flavour schemes of neutrino masses and mixing. Solutions of the solar neutrino 
problem are SMA (left panel) and LMA right panel) |0] 



on the vertical scale reflects their mass whereas the areas of differently marked parts are 
equal to |?7aiP, i-e. to the weights with which the ath flavor eigenstate is present in the 
ith mass eigenstate In the 3-flavour framework, the present-day data allow three main 
possibilities corresponding to the three main solutions of the solar neutrino problem. They 
are shown in fig. ^ and in fig. ^ (left panel). The schemes have different contributions of 
z/e, Ufj, and u-r in the two low- lying mass eigenstates. The weight of z/g in the mass eigenstate 
z/3 which is separated from ui and 1^2 by a large mass gap is |?7e3p; it is either small or zero. 
The weights of z/^ and z/,- within each mass eigenstate are approximately equal to each other 
as a consequence of 6*23 — 45° and |^7e3p = -C 1. This, in particular, means that the solar 
neutrinos oscillate into an almost equal superposition of and u^, independently of the 
type of the solution of the solar neutrino problem (SMA, LMA, LOW or VO). The scheme 



of the left panel of fig. ^ corresponds to single large mixing, whereas those of the right 
panel of fig. ^and left panel of fig. ^describe bi-large or bi-maximal (in the case Ue^ = 0) 
mixing. The 3-flavour schemes of figs |21| and ^ correspond to the so-called normal neutrino 
mass hierarchy, rrii ^ m2 ^ m-^. In addition, there are schemes (not shown in the figures) 
in which the position of the mass eigenstate z/3 and of the pair (z/i,z/2) are interchanged 
(inverted mass hierarchy, ms -C mi ~ ■^2)- The present-day data do not discriminate 
between the normal and inverted hierarchies; such a discrimination may become possible 
in future if the earth's matter effects in atmospheric or long baseline z/g ^ or z/g ^ Vt 
oscillations are observed. 

The normal and inverted neutrino mass hierarchies discussed above do not exhaust all 
possible schemes of 3-flavour neutrino mass schemes. Since neutrino oscillation experiments 
only allow the determination of neutrino mass squared differences and not of the masses 
themselves, it is also possible that neutrinos are quasi-degenerate in mass, and only their 



mass squared differences satisfy the hierarchies shown in figs. 21 and 22. Direct neutrino 



mass measurements allow the average neutrino mass as large as a few eV (provided that 
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Figure 22: Schemes of neutrino masses and mixing. Left panel: 3-flavour scheme with the VO 
solution of the solar neutrino problem. Right panel: SMA solution in the 4-flavour scheme [[fo|]. 



the 2j30i' constraint (|127| ) is satisfied). In that case neutrinos could constitute a significant 



fraction of the dark matter of the universe (hot dark matter). 

If the LSND result is correct, there should be at least four light neutrino species. Four- 
flavour neutrino schemes can be analyzed similarly to the 3-fiavour case. The data allow 
essentially two schemes. In each of these schemes there are two pairs of nearly degenerate 
mass eigenstates separated by large Am^^.^^. The mass splittings between the components 
of the quasi-degenerate pairs are AitLq and Am^^^. An example of a 4- flavour scheme is 
shown in fig. ^ (right panel). More detailed discussions of the 4-fiavour schemes can be 
found in . 

The experimental information on the neutrino masses and lepton mixing angles allows 
one to reconstruct the phenomenologically allowed forms of the neutrino mass matrix. This 
can be done by inverting the relation {mL)diag = diagimi, m2, m^) = U'^rriL U . The detailed 
structure of the neutrino mass matrix not yet known since some of the neutrino 

parameters are still rather uncertain. However, one can derive the zeroth-approximation 
structures (textures) of the neutrino mass matrices just using the already known gross 
features of neutrino spectrum and mixing angles. In particular, in the zeroth approximation 
one can assume mi = m2, sis = 0, 623 = 45°. Some of the neutrino mass matrix textures 
that can be obtained in this way are summarized in Table ^ taken from [|72[. Notice that 
these textures are defined only up to the trivial sign changes due to the rephasing of the 
neutrino fields. 

In realistic models, zeros in neutrino mass textures must be filled with small elements 
and in addition the large entries can be perturbed by small corrections. The small entries 
of the neutrino mass matrices should satisfy certain constraints; one of them follows from 
the requirement that for the LMA and SMA MSW solutions of the solar neutrino problem 
the lower-lying of the mass eigenstates ui and 1^2 have a larger u^. component. As was 
discussed in sec. ||, if this condition is not satisfied, the MSW resonance is only possible for 
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Table 2: Zeroth order form of the neutrino mass matrix for double and single maximal mixing. 
A - normal mass hierarchy, B - inverted hierarchy, C - quasi-degenerate neutrinos. From [^]. 
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antineutrinos and not for neutrinos. This requirement on the small entries of the neutrino 
mass matrices disfavours the inverted mass hierarchy and quasi-degenerate neutrinos and 
favours the normal mass hierarchy |]73 |. 

The neutrino mass matrix textures can provide us with a hint of the symmetries or 
dynamics underlying the theory of neutrino mass. With the forthcoming data from future 
neutrino experiments, it may eventually become possible to unravel the mechanism of the 
neutrino mass generation, which may hold the clue to the general fermion mass problem. 



We live in a very fascinating time for neutrino physics. It is very likely that in a few 
years from now new solar and long baseline neutrino experiments will bring us an important 
knowledge about neutrinos allowing to answer many questions and solving the solar neutrino 
problem. Neutrinos may also bring us new surprises, as they did many times in the past. 
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